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This  report  presents  major  accomplishments  and  significant  findings 
of  research  conducted  during  the  first  year  of  contract  to  the  FAA 
High  Altitude  Pollution  Program,  ending  June  10,  1976.  The  overall 

research  effort  is  concerned  with  numerical  modeling  of  the  atmos- 
pheric response  to  aerospace  utilization  of  the  stratosphere.  It 
encompasses  four  areas  of  study:  photochemical  kinetics;  coupled 

kinetics  and  transport;  radiative  transfer;  and  meteorological  anal 
ysls.  Additional  limited  studies  have  been  conducted  using  the 
Laboratory's  two-dimensional  climate  model  (ZAM2). 
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During  the  fiscal  year  ending  June  30,  1976,  luiwrence  Livermore 
Laboratory  (LLL)  participated  in  the  High  Altitude  Pollution  Program, 
which  was  instituted  by  the  U.S.  Department  of  Transportation's 
Federal  Aviation  Administration,  under  Reimbursable  Agreement 
DOT-TSC-76-1  with  the  Transportation  Systems  Center  in  Cambridge, 
Massachusetts.  This  report  describes  the  major  accomplishments  and 
significant  findings  of  the  LLL  work  performed  to  date.  The  overall 
research  effort  at  LLL,  which  is  primarily  concerned  with  numerical 
modeling,  has  been  divided  into  a number  of  subtasks  according  to 
research  area.  The  successful  accomplishment  of  these  subtasks  has 
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1.  INTRODUCTION 


The  Hisih  Altitude  Pollution  Program  (IIAPP)  was 
initiated  by  the  Federal  Aviation  Administration  for 
the  purpose  of  ensuring  that  stratospheric  aircraft 
engine  emissions  will  not  result  in  una  :.eptable  effects 
on  the  biosphere.  The  objective  of  H.APP  is  to 
quantitatively  determine  the  requirements  for  reduced 
cruise-altitude  e.\haust  emissions  and.  through 
international  and  interagency  cooperation,  to  ensure 
that,  if  necessary,  appropriate  regulatory  action  is 
taken.  HAPP  incorporates  many  related  activities, 
including  modeling,  field  measurements,  laboratory 
measurements,  engines  and  fuels  assessment, 
regulation,  and  monitoring. 

.As  an  active  participant  in  H.APP  during  the  fiscal 
year  ending  June  30.  1976.  Lawrence  Livermore 
Laboratory  (LLL)  has  undertaken  an  extensive  effort 
in  numerical  modeling  t>f  the  atmospheric  response  to 
aerospace  utilization  of  the  stratosphere.  A modeling 
effort  is  important  to  HAPP  for  two  primary  reasons. 
First,  modeling  constitutes  a means  of  relating  engine 
emissic'ns  to  changes  in  atmospheric  composition. 
Because  of  the  complex  physical  processes  and 
feedback  mechanisms  involved,  it  is  essentially 
impossible  to  assess  the  effects  of  emissions  without 
incorporating  these  processes  into  a model  and 
allowing  them  to  interact.  Second,  modeling 
contributes  to  the  interpretation  and  analysis  of 
experimental  results  and  helps  determine  the  relative 
priorities  for  carrying  out  this  work.  The  modeling 
effort  at  LLL  covers  four  major  areas  of  research: 
photochemical  kinetics,  coupled  kinetics  and  transport, 
radiative  transfer,  and  meteorological  analysis. 

Photochemical  kinetics  modeling  must  span 
gas-phase  reaction  processes  and  heterogeneous  (i.e., 
gas-particle)  reaction  processes,  constant  and  diurnally 
(and  seasonally)  varying  photcKhemical  reaction  rates, 
a wide  range  of  ambient  temperatures  and  pressures, 
neutral  and  ionized  species,  processes  with  both 
accurately  and  poorly  known  rates  of  reaction,  and 
species  of  both  acknowledged  and  qtiestionable 
importance  in  the  stratosphere.  Considerable  attention 
has  been  directed  toward  evaluation  of  the  sensitivity 
of  reaction  mechanisms  to  deficiencies  in  our 
knowledge  of  reaction  rates,  quanttim  yields,  and 
reaction  mechanisms. 

The  transport-kinetics  models  are  intended  for 
time-dependent  calculations  of  the  response  of  the 
stratosphere  to  injections  of  trace  species  by  aerospace 
vehicles.  These  one-  and  two-dimensional  models 
incorporate  interactive  radiative  and  chemical  processes 
coupled  with  transport,  which  is  described  by 
prescribed  mean  winds  and  eddy  diffusion  coefficients 


as  appropriate.  These  models  are  designed  to  facilitate 
iwrturbation  and  sensitivity  studies. 

Radiative  transfer  modeling  includes  ( I ) application 
of  highly  detailed  models  to  assessment  of  the  effect 
on  the  radiation  budget  of  perturbations  to  the 
stratospheric  composition,  and  (2)  development  of 
simplified  radiative  transfer  models  for  inclusion  in  the 
transport -kinetics  models.  Considerable  attention  has 
been  directed  toward  including  interactive  radiative 
processes  in  the  transport-kinetics  models  in  order  to 
assess  possible  feedback  mechanisms.  ' 

Meteorological  analysis  provides  guidance  and 
support  for  the  modeling  efforts  through  in-depth 
studies  of  atmospheric  processes  and  phenomena.  A 
primary  responsibility  of  this  research  is  to  develop 
test  situations  to  validate  various  aspects  of  the 
numerical  models. 

In  addition  to  these  major  research  areas,  special 
limited  studies  have  been  conducted  using  LLL's 
two-dimensional  climate  model  (ZAM2).  Climate 
model  development  at  LLL  is  presently  being 
supported  by  a related  project  involving  assessment  of 
the  clinratic  effects  of  atmospheric  nuclear  explosions. 
This  project,  which  is  funded  by  the  Division  of 
Military  Application  of  the  Fnergy  Research  and 
Development  Administration,  and  HAPP  both  take 
advantage  of  numerical  modeling  capabilities  and 
scientific  expertise  developed  at  LLL  during  the 
Climatic  Impact  Assessment  Program  (CLAP)  of  the 
Department  of  Transportation.  Although  the  same 
numerical  models  are  used  for  both  studies,  the 
applications  of  the  models  are  quite  different.  For 
example,  the  study  of  the  climatic  effect  of  nuclear 
explosions  is  primarily  concerned  with  the 
time-dependent  response  of  the  atmosphere  to  pulse 
injections  of  NO^.,  whereas  the  HAPP  study  is 
concerned  wiih  steady-state  injections  of  NO^.  Tasks 
relating  to  model  development  or  refinement  benefit 
both  projects,  and  consequently  such  tasks  have  been 
divided  between  the  two  projects.  In  this  way  each 
project's  sponsor  benefits  from  the  participation  of  the 
other  sponsor. 

These  two  projects  together  constitute  LLL's 
stratospheric  research  program.  The  various  research 
activities  included  in  the  program  and  their 
interrelationships  arc  diagrammed  in  Fig,  I . Fxcept  for 
the  two-dimensional  transport-kinetics  model,  we  are 
beyond  the  model  development  stage.  Because  of 
rapidly  changing  chemistry  and  rate  data,  most  of  the 
effort  in  the  transport-kinetics  area  during  the  past 
year  has  been  devoted  to  model  improvement  and 
sensitivity  studies.  Research  in  the  radiation  transfer 
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area  has  conccntraied  on  the  development  of 
subroutines  for  inclusion  in  the  transport-kinetics  and 
ciimate  models.  Several  sensitivity  and  validation 
studies  were  conducted  using  the  two-dimensional 
climate  model  (not  sponsored  by  IIAPP).  and  the 
results  from  these  studies  agreed  well  with  other 
models.  Meteorological  analysis  has  consisted  primarily 
of  comparing  tnodel  results  with  observations  when 
(X)ssible  and  developing  validation  tests  for  the  models. 

A long-range  objective  of  the  research  program  is 
the  development  of  a capability  for  assessing  the 
|x)ssible  climatic  effects  of  various  atmospheric 
ixrturbations.  Currently,  our  assessments  of  perturbed 
sivcies  concentrations  consider  a large  number  of 
photochemical  reactions  and  are  capable  of  treating 
radiative  and  diurnal  processes  in  detail,  and  they 
assume  that  the  effects  of  transport  on  atmospheric 
composition  can  be  represented  parametrically.  Our 
climatic  assessments  treat  many  interrelated 


thermodynamic  and  hydrologic  processes,  and  the 
ix'rturbed  chemical  composition  is  treated  as  a given 
t|uantity.  Ultimately,  our  assessments  will  include 
interaction  between  the  chemistry  and  the  dynamics. 

The  sensitivity  studies  conducted  during  the  past 
year  were  intended  to  increase  our  understanding  of 
the  characteristics  of  the  numerical  models  and  the 
uncertainties  associated  with  their  results.  The 
ix'rturbations  (e.g.,  injections)  used  in  these 

studies  were  chosen  so  as  to  cause  a significant  model 
response  rather  than  on  the  basis  of  an  expected  .SST 
fleet  size.  The  reader  is  cautioned,  therefore,  not  to 
interpret  the  results  as  predictions  of  the  effects  of 
aerospace  operations,  but  rather  as  indications  of  the 
model  sensitivity.  Predictions  of  the  effects  of 
aerospace  operations  on  atmospheric  composition  must 
include  the  simultaneous  injection  of  several  species 
(e.g..  NO.,^.  H-,0.  and  SO-,).  During  the  next  yeai 
simultaneous  injections  will  be  considered. 
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Fig.  1.  Flow  chart  showing  the  approach  I.l.L  ha.s  taken  in  its  stratospheric  research  program. 


2 


2.  RESEARCH  ACTIVITIES 


1 


Because  of  the  complex  nature  of  the  work,  most 
of  the  research  tasks  represent  the  combined  effort  ol 
several  scientists  ratlier  than  individual  contributions. 
It  is  not  possible  to  separate  the  tasks  by  research  area 
because  of  the  team  approach.  Major  tasks  that  have 
been  completed  are  reported  here  in  considerable 


detail.  Work  still  in  progress  and  several  shorter  tasks 
that  have  been  finished  are  summarized  at  the  end  of 
this  section.  Those  contributions  which  we  feel  are 
particularly  significant  have  been,  or  are  in  the  process 
of  being,  submitted  for  publication  in  technical 
journals. 


2.1  NO^  CATALYTIC  OZONE  DESTRUCTION:  SENSITIVITY  TO  RATE  COEFFICIENTS 


In  recent  years  it  has  been  suggested  that  several 
of  man's  activities  could  result  in  a reduction  of 
stratospheric  ozone.  This  would  increase  the  UV-B* 
l1u.\  at  the  earth's  surface,  and  presumably  result  in 
an  increase  in  the  incidence  of  skin  cancer.  It  has  also 
been  suggested  that  changes  in  stratospheric  ozone 
might  alter  the  climate  (Grobecker  et  al.  I974). 
.Among  others,  from  stratospheric  aircraft  tlights 
j (Johnston  1971,  Crutzen  1972).  NO,,  from  nuclear 

I tests  (Foley  and  Ruderman  1973).  and  nuclear  warfare 

' (MacCracken  and  Ctiang  1975.  National  Academy  of 

, Sciences  1975a)  have  been  proposed  as  capable  of 

reducing  stratospheric  ozone.  In  order  to  estimate  the 
I magnitudes  of  these  hypothesized  effects,  considerable 

' effort  has  been  devoted  to  the  construction  of 

^ numerical  models  of  the  stratosphere  (Grobecker  et  al. 

1974.  National  Academy  of  Sciences  1975b. 
Interdepartmental  Committee  for  Atmospheric 
, ' ^ Sciences  1975)  and  acquisition  of  data  to  be  used  for 

model  input  and  to  verify  model  results.  The  major 
reasons  for  current  credence  in  the  nuxlels  have  been 
•j  that  they  have  used  measured  inputs  (so  that  the 

components  of  the  mode!  are  presumed  to  be 
•!  separately  tested)  and  that  the  models  generate 

concentrations  of  stratospheric  species  that  fall  within 
the  sometimes  rather  wide  range  of  measured  values. 
• While  the  above  considerations  are  necessary  for 

J model  credibility,  they  are  not  conclusive  for  several 

; reasons.  First,  the  models  are  known  to  be  incomplete 

j in  various  ways.  All  use  more  or  less  severe 

( approximations  in  treating  photon  Oux  densities.  Ail 

° 1-0  and  2-D  (one  dimensional  and  two-dimensional) 

« models  use  parameterized  transport.  All  the  models  use 

an  incomplete  representation  of  the  chemistry  (for 
I example,  ionic  processes  and  heterogeneous  processes 

■ are  omitted  or  very  crudely  represented  in  existing 

1 models).  .Second,  comparison  of  computed  and 

I observed  concentrations  is  subjective  since  most 
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parameters  reveal  a substantial  range  of  measured 
values  and  hence  some  selection  must  be  made  of  the 
actual  data  to  be  compared.  Finally,  each  of  the  input 
parameters  has  some  uncertainty  associated  with  it 
(e.g..  in  the  case  of  the  chemical  rate  data  the 
uncertainties  that  have  been  estimated  for  many  of  the 
processes  range  from  a few  percent  up  to  a factor  of 
10  (Uampson  and  Garvin  1975)). 

The  first  two  problems  mentioned  are  in  a sense 
philosophical  and  may  not  be  addressable  in  the 
framework  of  existing  "mean  state"  models  with 
admittedly  restricted  degrees  of  freedom.  The 
estimated  effects  of  processes  known  to  be  omitted 
are  a few  percent  or  less.  Although  ionic  processes  have 
recently  been  suggested  to  be  important  (Ruderman 
et  al.  1976),  this  argument  is  speculative.  However,  the 
third  source  of  uncertainty  can  he  examined: 

In  principle  there  are  several  techniques  available  to 
study  quantitatively  the  propagation  of  error  in  a 
complex  model  (Dickinson  and  Gelinas  1975,  C'ukier 
et  al.  1973.  Schaibly  and  Schuler  1973.  and  .Atherton 
et  al.  1975).  Linfortunately . at  present  these  techniques 
are  prohibitively  expensive  and  are  difficult  to  apply. 
Nonetheless,  it  is  feasible  to  perform  experiments 
which  examine  model  sensitivity  to  perturbations  in 
individual  rate  coefficients  or  ensembles  of  rate 
ctrefficients  chosen  with  the  intent  of  maximizing  or 
minimizing  the  model  response  to  an  incremental 
change  in  some  other  parameter.  This  would  then 
provide  approximate  bounds  for  the  model's  sensitivity 
to  rate  coefficient  uncertainties.  In  this  study  we  adjust 
the  rate  constants  involved  in  a tmxlel  of  the  HO.,^. 
NOj.,  and  reactions  in  the  stratosphere  and  then 
exatnine  the  effects  on  the  response  of  the  model  to 
an  NOj^  injection.  This  technique  is  useful  in 
determining  extreme  ranges  of  possible  effects,  and  in 
identifying  highly  sensitive  parameters  ot  groups  of 
parameters.  It  is  recognized  that  this  is  not  very  useful 
in  estimating  the  probability  of  errors  of  any  given 
magnitude.  But  it  could  give  qualitative  lower-limit 
estimates  of  the  expectation  of  an  el'fect.  FDr  even 
this  to  be  done,  it  is  necessary  to  estimate  the 
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ilisinbiilion  of  expected  error  in  lire  individual  rate 
eoelTieients  and  the  sensitivity  of  ilic  nurdel  to  subsets 
ol  the  rate  eoelTieients. 

The  i|iiantitative  nteaiting  of  the  nireertainties 
i|iioted  by  llaittp.soir  and  (iarvin  (ld75)  aird  used  iir 
tlris  stndv  is  not  easily  speeified.  The  ipiotcd 
nneei tainties  were  inteirded  to  defiire  a range  that 
would  be  ex|X'eted  by  a panel  of  experieneed  ehentieal 
kinetieists  to  eontain  the  eorreet  valne  (llatitpson 
We  have  exaittined  data  pirblished  after  tire 
evalnations  given  by  Haitrpson  and  (iarviir  for  six 
leaetions  of  partieular  ittrpoitanee  to  our  nrodel  and 
sotite  0(*D)  reactions.*  Wc  find  that  foi  krwer 
stratospheric  conditions  the  mean  divergence  between 
new  data  and  the  value  cited  by  llampson  and  Garvin 
was  about  1.2  times  the  cited  uncertainty,  but  that 
the  error  distribution  was  very  different  from  a normal 
distribution. 

If.  in  spite  of  tire  above  finditrg.  one  treats  the 
quoted  uncertainties  as  estiitiates  of  the  expected 
errors  (i.e.,  I a),  then  the  expectation  that  all  of  the 
varied  rate  coefficients  will  depart  at  least  as  far  from 
the  "best  " values  as  those  used  in  a test  model,  and 
in  the  same  sense,  will  be  roughly  (0.16)"  where  n 
is  the  number  of  adjusted  rate  coefficients.  However. 


•The  reactiuns  examined  were  Iv6.  KIO,  Kll.  kl2.  K14, 
K18,  K19.  K20.  K26.  and  the  reaetions  of  ()('!))  wiili  Ni, 
Ol.  and  COj.  OnK  relative  rates  were  enmpared  Inr  the  C)( ' I)) 
reaetions  with  speeies  other  than  fOi.  which  were  the  stuiree 
ot  absolute  rate  data  in  llampson  and  llarvin 


because  the  tnodel  is  insensitive  to  many  of  the  rate 
coefficients  and  ver\  sensitive  to  a few.  this 
methodology  piovides  only  a lower  limit  to  the 
expectation  of  a model  with  a sensitivity  at  least  as 
divergent  from  the  "best"  model  as  the  model  tested. 
If  many  of  the  parameters  varied  are  insensitive,  this 
lower  linrit  would  be  so  extreme  as  to  be  meaningless. 

The  basic  model  used  for  these  experiments  was  the 
LLL  one-dimensional  model  used  in  the  (TAl’ 
dticumentation  ((jrobeckei  ct  al.  l‘>74.  Grobecker 
1675).  Its  structure  has  been  described  elsewhere 
(Chang  et  al.  1674).  This  model  includes  the  chemical 
rcactitms  given  in  Table  1.  It  tieats  the  lolUiwing 
s)x;cies  dynamically  (by  parameteii/ed  vertical 
transport)  as  well  as  photochemicalK : O.  (),.  NO. 
NO,.  UNO,.  110.  HO,.  11,0,.  N,0.  Steady-state 
approximations  are  made  for  ()(*D).  N.  and  11.  and 
the  concentration  profiles  of  11,0  and  CH^  are 
prescribed.*  The  various  species  concentrations  are 
.solved  numerically  at  44  levels  extending  to  an  altitude 
of  35  km.  In  this  study  both  the  carbon  cycle  (except 
insofar  as  methane  is  a source  of  110^  ) and  the  halogen 
cycles  were  omitted.  The  response  of  the  model  to 
the  choice  of  eddy  diffusion  profile  w as  examined  bv 
Chang  (1674).  For  this  study  we  used  the  LLL 
dilfusion  profile  derived  in  the  earlier  study  (Chang 
et  al.  1674).  (See  Fig.  2.)  Nine  sets  of  rate  coefficients 
were  examined: 


•This  version  ot  the  basic  model  is  elioscn  speeifiealh  to 
provide  continuity  with  the  CT.AP  Report. 


Table  1.  Reactions  and  rate  parameters  for  the  various  models  used. 


I’hol'rlysis  rates  al  the  indicated 

Keaetiiin  altitudes  for  Model  .Aj 


Formula  Number  55  km  30  km  IS  km 


(>2  + hr  -* 

20 

.11 

7.76 

X 

10- "> 

1.90 

X 

10" 

8.42  X 

10'*’ 

(>3  -r  h.  • 

0 + Oj 

.12 

2.80 

X 

lO"* 

2.41 

X 

in'* 

2.07  X 

10-^ 

O3  + hi’  • 

0(’0)  + O2 

.13 

4.79 

X 

10'^ 

5.85 

X 

10'^ 

7.90  X 

10'^ 

NOj  + he 

• NO  + 0 

.14 

5.03 

X 

10'^ 

4.85 

X 

10'^ 

4.74  X 

10'^ 

NjO  • hr 

N^  + O(’O) 

.15 

4.14 

X 

10-^ 

2.89 

X 

in** 

1.44  X 

lo'-'’ 

NO  + he  - 

N + 0 

.16 

I..35 

X 

Hf^ 

8.90 

X 

10  "> 

7.30  X 

10-’’ 

IINO3  + hi 

- OK  + NO3 

.17 

6.38 

X 

10-5 

4.78 

X 

10*’ 

3.21  X 

10’ 

lljOj  + hi 

■ 2011 

IS 

8 25 

X 

10-' 

6.39 

X 

lO*’ 

2.22  X 

10'’ 

IIO3  + hi 

- Oil  + 0 

.14 

1.98 

X 

itr‘ 

.3.40 

X 

10*’ 

1.70  X 

10"' 

NO3  + hi 

■ NO  + (>2 

III), 

(1  ilht’r 

.II0|  ot.l|02  is 

assu  met 

In  be  Ihc 

only  fate  of 

NO  3 ♦ hr 

■ NOj ♦ 0 

'l"2 

N03 

as 

IS  inJii  atcil  hv 

imulol 

siibscrip*.) 
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1 i>ns(.irit  \ m K \ o\p(<  I ) in  unils  ol  cm  /s  (cm  /\  for  thrcc-lindy  processes), 

onsi.int  (.  in  K ,\  exp(l.  < I ). 


Model  X'  dilTers  from  \ tmls  in  (h;il  KI9  is  sel  t»»  2 X 10*  . (he  vjiue  used  in  nmst  ol  (he  C'lAI'  expenmenls 

'^K2()  IS  pressure  de|H*ndent.  In  our  model  K20  K|M|/(I  12  X l()'^  + |X1|)  where  K is  \;ilue  tabled  above.  Ibis  is  a fit  to  the 

altitude  tiependen I eHpressi«ui  uiven  bv  llainpstm  and  (iarvin  (l'^^?). 


r 


Fig.  2.  profiles  that  have  been  used  by  various 

investigators.  The  Oiang  (1974)  profile  was  used  to  examine 
the  nine  sets  of  rate  coefficients  described  in  the  text. 


1.  Model  A I used  the  rate  eoelTicients  given  In 
Hainpson  and  (iarvin  (1975)  with  llie  rate  eoelTicient 
lor  the  reaction  HO  + HO,  • 11,0  + O,  (KI9)  set 
to  2 X 10"'  ' etna  's.  Tliis  is  the  low  end  of  the  rattge 
i|uoted.  and  the  value  currently  preferred  by  us 
(Kaufman  1975.  I.loyd  1974).  Branch  I was  assumed 
for  J 1 0. 

2.  Model  A,  differed  from  model  .Aj  in  that 
branch  2 was  assumed  tor  JiO. 

,5.  Model  A 'j  differed  from  model  Aj  it)  that  K19 
was  set  to  the  high  end  of  the  range,  i.e..  2 X lo  ''*. 
the  value  used  for  most  I'f  the  ( lAP  computations. 
Branch  I was  assumed  for  JIO. 

4.  Model  A differed  from  model  A ^ by 

assuming  branch  2 for  JIO. 

5.  Model  B assumed  all  tlie  rate  coefficients  were 
in  error  by  the  estimated  liticeftaitity  indicated  by 
Hampson  and  (iarvm  (1975)*  and  were  varied  within 
this  range  with  the  intent  of  maximizing  the 
destruction  of  ()j  by  a given  atnouiii  of  NO.^. 
minimi/ing  the  destruction  of  (),  by  ambient  and 
IIO^.  and  miitimi/ing  the  ambient  concentrations  of 
\()^  and  IIO^.  Branch  I was  assumed  lor  JIO. 


'When  a r.itc  cnnsl.int  wav  reported  to  be  known  wjih  higlur 
accuracy,  but  vmly  at  room  tcmpcralurc.  and  tlic  roenn* 
temperature  rate  cimsi.int  w av  lew  rlun  10*7  cm-^'v.  .i  '3-told 
uncertainty  wav  avvumed. 


I 

L 


6.  Model  ( | varied  all  rate  coeKicients  withiti  the 
same  range  of  unceitainty  with  the  intent  of 
minitnizitig  the  destiuctioii  of  by  a given  c|uantiiy 
o!  NO^.  maximizing  destruction  by  ambient  ()^ 
and  IK)^.  and  tnaximizing  the  ambient  concentrations 
of  NO  and  HO  . Bianch  I was  assumed  for  JIO. 

X \ 

7.  Model  {’,  dilTered  fiom  model  ( j by  assuming 
branch  2 for  jFo. 

b.  In  an  attempt  to  identily  the  leactions  most 
responsible  lor  the  chatige  m model  behavior,  a model 
(' I . analogous  to  model  Cp  was  developed  in  which 
only  a limited  set  of  rate  coefficients  was  adjusted. 
Brandi  1 was  assumed  for  JIO. 

9.  Model  ('  , ditTered  from  model  ('  | by 
assuming  branch  2 lor  JIO. 

It  should  be  noted  that  the  uncertainties  in  the 
photodissociatioii  rates,  wliich  are  probably  even  larger 
than  the  uncertainties  in  thermal  rate  coefficients,  were 
not  treated  here.  Several  approximations  are  involved 
in  the  calculation  of  photodissociation  rates,  and  these 
approximations  are  the  subject  of  concurrent  research. 

The  following  computations  were  carried  out: 

1.  An  ambient  cotidition  was  cotiiputed  for  each 
of  the  models. 

2.  The  condition  after  .H)0  years  of  NO,  injection 
ai  2000  tiiolecules/cttr ‘s  throughoui  a I -knt-thiciv 
layer  at  17  km  altitude,  the  average  cruise  altitude  of 
the  Concorde  SST.  hereafter  designated  as  NO^-17.* 

.5.  All  models  were  used  to  compute  the  condition 
after  300  years  with  the  above  NO,,  injection  at 
20  km,  a possible  cruise  altitude  for  an  "advanced" 
SST  (NO,^-20). 

4.  Same  as  cases  2 and  3 w ith  injection  at  35  km 
(NO,^ -35). 

5.  Antbien)  conditions  (case  I ) were  recomputed 
(or  all  models  with  the  water  column  increased  by 
lO';.  roughly  the  magnitude  of  the  increase  in 
stratospheric  water  content  exjx’cted  from  the 
hypothetical  fleet  of  500  Boeing  SST's  (Grobecker 
1975)  (MX  11,0). 

6.  Case  2 was  repeated  with  the  water  column 
increased  by  107;  (1.1  X 11,0.  NOpl7), 

The  calculated  Oj  columns  atid  the  percetit  of  O, 
deviation  (AO,)  front  the  respective  anibietit 
conditions  are  given  in  Table  2.  The  effect  of  the 
choice  of  diffusion  profile  was  examined  for  cases  1 
and  2 and  the  results  are  summarized  in  Table  3.  The 
piofiles  studied  are  given  in  l ig.  2. 


'ITiiv  emivsion  rate  eorrevpondv  to  a very  larec  Coneortle 
nvx't,  mmibcrcii  in  thoiivamls  ot  planov.  ( or  eomparivon  we 
examined  a tew  eases  assuming  a leni'old  smaller  mieetion  rale 
and  foiind  llie  resultant  pcrUirbations  to  be  about  ninelold 
smaller,  suggesling  an  almost  linear  relationship  between 
inieeiiim  rate  and  resultant  perturbation. 
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Table  2.  Calculated  ambient  and  perturbed  column  densities  (10’®cm  and 
percent  deviation  from  the  respective  ambient  column  densities 

for  the  perturbation  cases  1-6  as  described  in  the  text.  See  text 
for  description  of  models  ■ ■ — 


Case 

(1  \P  llUHld 
(hi);h  value  of  K 19) 

*'Best”  model 
(lou  value  of  K 19) 

Mavimi/ed  re- 
sponse to  NO^ 

Miriiim/ed  o'sponse 
lo  NO 

Response  6)  NO^  mini 
mi/ed  on  four  rales 

'i 

'2 

A, 

B 

7 

C, 

‘ 1 

*^2 

JIO  yields 
NO  Oj 

J 10  yields 
Nt>2  + O 

J 10  yield^ 
NO  + (>2 

JIO  yields 
NOj  + O 

JIO  yields 
NO  + (>2 

TuTyidds 
NO  + (>2 

JIO  yields 
N(>2  + O 

J 10  yields 
NO  + (>2 

J 10  yields 
NO2  + If 

(.  Amhtcnr 

7. 765 

7.993 

7.422 

7.555 

7.577 

6.893 

6.997 

7.338 

7.415 

2.  NO  17  0 

7.4  2H 

7.712 

7.292 

7.467 

6.996 

7.013 

7.152 

7.383 

7.490 

4.347. 

3.52'i 

■1.757 

1.167 

■8.067 

+ 1.747 

+2.227 

+0.617 

+ 1.017 

3.  20  <)j 

7.005 

7.316 

7.035 

7.250 

6.370 

7.004 

7.182 

7.323 

7.467 

.sDj 

•9.79'/^ 

«.477- 

■5.217, 

4.037 

15.927 

+ 1.617 

+2.647 

-0.20", 

+0.70'/, 

4.  NO^  35  Oj 

6-426 

6.767 

6.377 

6.614 

5.674 

6.450 

6.663 

6.663 

6.837 

50, 

■17.24'^ 

15.347 

14.067 

12.467 

25.117 

■6.4  2'-r 

4.777 

•9.20'<'r 

■T.my/r 

5.  I.l  X 11,0  0, 

7.767 

7.992 

7.408 

7.539 

7.588 

6.857 

6.959 

7.312 

7.387 

+0.037r 

0.017 

0.147 

■0.217 

+0.157 

0.527 

0.547 

■0.357 

-0.387 

6.  I.l  X 11,0.  0, 

7.435 

7.715 

7.282 

7.454 

6.979 

6.9"8 

7.114 

7.357 

7.462 

NO^  17  50, 

•4.25"r 

-3.48'I 

1.887 

-1.347 

-7.89'^ 

+ 12.(3 

+ 1.677 

+0.267 

+0.637 

Tables.  Calculated  ambient  (case  1)  and  perturbed  (case  2)  column  densities  IW^^cm'^)  and  percent 
deviation  from  the  respective  ambient  column  densities  for  four  eddy  diffusion  coefficient 

profiles  as  described  in  the  text.  See  text  for  description  of  models  


I ddy  diffusion 
aH.’fficient  profilo 


Model  t'2 


Model  A 


Model  A 


Model  B 


Model  C 


Altitude  - km  Altitude 


0^  data: 


30 1—  • Model  A 

! □ Model  B 

o Model  C 
1 Krueger  (1973) 


20 


lOh 


10 


10 


10 


J1 


13 


0^  concentration  - molecules/cm'' 


|.  3.  (alciilated  equilibrium  concentration  profiles  for  compared  with  observations  from  Krueger  (1973).  Note  that  for 
igs.  3'7  all  surface  concentrations  are  input,  not  calculated. 


60  ( 


" I T 


50 


(Use  upper  scale) 
(Use  lower  scale) 


40  h 


NO^  data: 


30  h 


• Model  A 
o Model  B 
o Model  C 


20 


Ackerman  (1975) 

(data  summary)  _#■ 


10 


Ol 


I 1 


•-ad 


ci-n^.1 


10  10"  10"  10'"  10 
NO2  concentration  — molecules/cm 

Fig.  4.  (alculated  equilibrium  concentration  profiles  IV»r  Ml,  compared  with  data  •aimmary  (d  Ackerman  (197^). 
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1 ho  ooin|Hiled  ainbiciil  oimooiilralion  profiles  aie 
ooiiipaiod  with  observations  in  l ies.  3-7  for 
models  A|.  If.  and  ( (l-or  most  speeies  model  A' 
closely  resembled  model  H and  model  C'  resembled 
model  t'.  Differences  due  to  the  JIO  branch  were 
baiely  detectable  on  tbese  semiloj;  plots,)  Only  profiles 
tor  s|X'cies  for  which  we  could  find  measurements  are 
ijiven  here, 

,As  can  be  seen,  the  computed  ambient  states  tor 
the  species  that  have  been  measured  in  the  stratosphere 
are  rather  similar  for  all  three  models  (labeled  simply 
as  A.  B.  and  ('  in  the  figures).  The  concentrations 
plotted  represent  dmrnally  averaged  conditions  as 
parametei  i/ed  by  a "half  sun"  condition  (the 
photoKsis  rates  used  are  one-half  of  the  noontime  rate 
as  calculated  by  the  model),  Oitly  for  the  OH 
measmements  of  ,\nderson  (Ih^p)  ^;J|l  crrmputed 
concentration  protiles  be  said  to  disagree  with 
measurement  by  mote  thatt  tbe  spread  in  recent 
measurenrents,  I rom  l ig.  ^ rtiod’l  If  would  be  the 


most  likely  to  be  judged  inadeiittate  However,  in  order 
to  be  more  detimtive.  mie  must  carry  out  detarled 
sensrtivity  analyses  insolving  ditrrnal  rtrodels  on  the 
details  of  HO^  chemrstry  sittce  OH  vat  rations  are  Inghly 
tirne-deixMident 

Based  on  presetrt  krrtetics  data,  ntodels  B arrd  ( ate 
highly  improbable.  The  fact  that  they  gave  contpttted 
profiles  agreeing  as  well  as  rttcxlel  ,\  with  present 
observatirrns  irrdicates  that  such  agreement  is 
insufficiertt  to  validate  a nuKlel,  This  seems  to  sttggest 
that  lor  observational  data  to  be  ttsed  to  dtserirninate 
atttong  models  the  concentratiotrs  must  he  kttown 
within  a factor  of  two  at  the  very  least.  If.  as  seems 
probable,  rttany  trace  species  vary  strhstantially  in  time 
and  space,  the  above  statement  pertains  to  the 
aeeitracy  with  which  the  correentration  distribution 
rnirst  be  knowtt.  a much  more  e.xpensive  undertaking 
than  a single  rneasurenrent. 

Berttrrhatirrns  in  the  O^  crriumn  computed  for  a 
17-km  injectiorr  of  NO^  (case  3)  are  plotted  ott  a 


10" 


10- 


10 


10 


NO  concentration  - mo1ecu1es/cm- 


Fig.  5.  Call.  It  la  led  eqiiilihriiim  eoneeniralion  prol'iles  for  NO  compared  with  data  summary  of  Ackerman  (1975). 


liiKMi  scale  III  I le.  As  can  lie  seen,  llie  liopos|i|ieiic 
[X'l till bainms.  laiitely  inniieiiced  In  the  ( ll_j  and  N()^ 
sinov;  reactions  (Cliaineides  and  Walkei  Johnston 

and  Qiiitevis  l')74).  aie  i|inte  small  coiiipaied  to  the 
stiatospheric  pei tin hatioiis.  This  is  in  part  due  to  out 
use  ot  li\ed  coiiceiiliations  as  a siiilace  hoiindais 
condition.  We  have  e\[ieiiineiited  with  the  lower 
houndaiA  condition  stilTicientIv  to  deinonstiate  that 
this  t\  pe  ot  surtace  boinulaiy  condition  has  a nejtligible 
elTeci  on  computed  stratos|ihertc  profiles.  althoui;h  the 
trojuisphertc  sensitivities  aie  sijtnilicantlv  ailected  by 
the  boundary  conditions. 

The  opinion  has  been  expressed  that  model 
piedictions  ol  (),  destiiictioii  loi  a eiven  N()^  miection 
are  accurate  to  within  a tactor  ot  2 or  .1  ((iiobeckei 
et  al.  l'-)‘'4).  ,As  can  be  seen  m Table  2.  the 
ix'rturbations  iti  the  ()..  coluitm  calculated  liotn  the 
base  case  (model  A I and  the  most  sensittve  case 
(tnodel  B)  dtl'fered  by  a factor  of  roughly  5 for  a 
17-kni  injection.  .J  for  a 20-km  injection,  and  only  I .(1 
for  a .s5-kni  injection.  Tims,  m the  case  of  ma.\imi/ing 
the  effect  of  an  NO^  injection  the  uncertainty  is 
compaiable  to  that  expected,  atul  is  not  severe. 
Tuithei.  more  than  half  the  cotnputed  difference 
resulted  from  the  variation  of  K Id  (compare  models  A 


and  A’Ioi  case  2).  Ilowevei.  in  (he  case  ol  niinimi/ing 
O.J  destruction  (model  ( ).  an  N()^  injection  at  17  km 
resulted  in  an  iiurcasc  m ()^  ei|ual  to  about  2.‘i''i  ol 
the  deciease  in  me  base  case.  Toi  a .J.'i-km  inieclion 
ol  NO^  all  of  the  models  predicted  an  Oj  deciease. 
the  model  C decrease  being  about  40'/f  as  laige  as  thai 
of  model  A.  Thus,  the  range  ol  piedicied  ellecis  is 
small  for  a high  altitude  injection,  and  the  unceilainiy 
is  modest  foi  maxinii/mg  the  effect  ol  an  \().^ 
injection;  however,  the  uncertainly  is  very  substantial 
(encompassing  a reversal  ol  the  el  led)  in  the  case  ol 
minimi/ing  the  effect  ol  a given  NO^  inieclion  at  low 
altitudes. 

If  it  were  necessary  to  vary  all  .1.1  ihernial  late 
coefficients  to  the  limit  of  their  unceilainty  in  oidei 
to  obtain  this  reversal,  the  expectation  of  its  being 
realized  would  be  almost  nonexistent.  However,  we 
examined  the  results  of  the  model  ('compulation  and 
were  able  to  identify  several  key  parameters  that  aie 
responsible  for  the  reversal  of  the  calculated  efiecl. 
We  found  in  model  ('  that  al  low  altitudes  ()., 
is  regulated  primarily  by  the  110^  cycle  and  that  the 
most  inipoilant  parameters  were  (hose  regulating  the 
110^  concentrations,  the  reactions  ol  110^  with  odd 
oxygen,  and  the  reactions 


-s; 
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HNO^  concentration  - molecules/cm" 


Fig.  6.  (alciilatcd  egiiilibriiim  eoiKTiilraliiin  profiles  for  IINO^  compared  with  data  summary  of  Aekerman  (iy75). 
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Fig.  7.  Calculated  equilibrium  concentration  profiles  for  Oil  compared  with  observations  by  Anderson  (1976)  and  Davis  (1975). 
Most  of  Andetstjn’s  data  were  talten  at  high  /enitli  angles,  whereas  Davis's  data  were  taken  at  low  zenith  angles.  Thus,  the 
“half  sun"  correction  should  not  be  applied  in  making  comparisons  with  (he  Anderson  data. 


IK),  + NO- NO,  + 110  (K:6) 

ami 

M 

HO  + NO,  -»  IINOj.  (K:0) 

These  reduce  the  elTicicncy  of  the  NO^  ami  IIO^ 
cycles  in  the  lower  stratosphere  where  the  NO^ 
rlestniction  cycle  is  already  inefficient  because  NO, 
is  usually  phololy/.cd  there  resulting  in  the 
regeneration  of  odd  oxygen. 

To  improve  our  understanding  we  have  sought  to 
sepaiate  the  odrl  oxygen  destruction  loles  of  the 
Chapman.  N()^ , and  HO^  cycles.  Since  the  cycles  are 
known  to  he  coupled  (e.g..  the  NO^  and  IIO^  cycles 


are  coupled  through  K20.  K21.  and  K2(i).  it  is  unclear 
which  rpianiiiies  would  be  most  illuminating.  As  a 
beginning  we  have  defined  the  net  odd  oxygen 
destruction  rate  (D)  of  each  of  the  cycles  as  follows: 

D(Chapman)  = 2K2|0|  [Ojl. 

(O,  photolysis  is  omitted  because  we  wish  to 
compare  desi ruction  terms.) 

IXNO^i  k.l  |NOHO,|«  M lNO,l|0|»  s:5  |N01|O)|M|->  k."l  INoplO,) 

J4  |Nop  Is  |NO|  kl  |N||Oj|  K»  |N|  |N()|  K.'4  |Nl|NOj| 

* JU)  (NO,  |(ii  Branch  I (' asMitncd) 
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(All  reactions  involving  jN]  are  trivial  below 
50  km.  as  is  also  NO  photiilysis;  JI0|N()^|  is 
treated  as  e<,|ual  to  the  rate  of  NO,  formation.) 

rvrroj  Kr:  |o,|iotrit  K!.i  |o|lori|t  Kr4  io,ntK>,|*  Krsioi  |tro,| 

♦ rsr7  |0|||rr|t  Kro  |o'i)| |rr.o|  rc:7'|orr|-'  |rro,| 


(The  last  two  terms  are  of  little  significance.) 

Vertical  profiles  of  these  quantities  evaluated  under 
ambient  or  equilibrium  conditions  of  models  A,.  B. 
and  C^  are  compared  in  Fig.  0.  The  largest 
differences  between  tlic  models  appear  in  D(HO^), 
in  D(C'hapman)  above  40  km,  and  in  D(NO,^) 


below  27  and  above  45  km.  Table  4 was  prepared  to 
reflect  these  differences  qirantilatively.  It  gives  for  each 
model  the  relative  traction  of  the  integral  (over  the 
altitude  range  indicated)  of  the  odd  o.wgen  destruction 
rale  ascribable  to  each  of  the  three  cycles. 

The  model  most  comparable  to  the  models  used  in 
ClAF  is  A '] . The  instantaneous  rates  widely  quitted 
from  Johnston  (l‘)75)  are  included  for  comparison. 
The  differences  in  results  retlect  differences  m 
methodology  and  definitions  as  well  as  differences  in 
rate  coefficients  and  species  concentration  profiles. 
Our  results  were  obtained  from  a 1-D  time-dependent 
model  with  coupled  kinetics  and  transport  e.xtending 
from  0 to  55  km  using  photolysis  rales  averaged  over 
the  diurnal  cycle  by  taking  halt  the  rates  calculated 
for  a 45°  /enith  angle.  Johnston's  values  were  obtained 


AO^  - 10^^  molecules/cm^ 

Fig.  8.  Calculat'd  pcriurbalions  in  (Ik*  O ^ pronic  IV>r  a I T-km  injcctUtn  of  NO^  (case  2>  plotted  as  departures  from  the  ambient 
profiles  (case  It  of  the  models  indicated. 
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Fig.  9.  Calculated  equilibrium  profiles  of  odd  oxygen 
destruction  rates  as  defined  in  text  for  ntvsdeis  \,  ii,  and 
C.  Where  D(.NO^)  forms  cusps  at  the  left  border,  it  changes 
sign.  Regions  of  net  ozone  production  by  NO^  are  shaded. 


from  a itoiulynumic'  iiukIcI  integrated  hetweeii 
"Impopause  ' and  4.‘>  km.  but  using  photolysis  rales 
lor  siK'eilie  altitudes,  laliludes.  and  times  of  day.  Little 
t|uantitative  significance  should  be  attached  to  these 
integrated  destruction  rales  since  they  do  not 
accurately  rellect  the  relative  importance  of  the  various 
processes  in  determining  total  O^.  In  an  effort  tr> 
estimate  the  relative  importance  of  these  processes  in 
determining  total  column  densities,  we  have  weighted 
the  various  local  destruction  rates  by  the  product  of 
the  local  Oj  concentration  and  ihe  local  lifetime*  for 
each  level  and  integrated  the  results  from  I I to  km. 
These  weighted  loss  rates  are  insensitive  to  the  choice 
of  the  integration  limits  so  long  as  the  top  is  above 
45  km  and  the  bottom  is  below  12  km. 

In  an  effort  to  identify  the  most  sensitive  reactions 
we  generated  additional  models  by  varying  smaller 
subsets  of  the  33  K‘s  listed  in  Table  1.  We  found  that 
varying  as  few  as  three  rates.  K14.  KIb.  and  K26.  was 
sufficient  for  an  NO^  injecliiin  at  17  km  to  produce 
an  increase  rather  than  a decrease  in  the  column. 
Model  C uses  the  values  of  Ihe  three  K's  cited  above 
plus  K20  from  model  ('  and  the  rest  from  model  A. 
As  seen  in  Table  2 this  is  sufficient  for  a 17-ktn 
itijeciioti  of  \'0^  to  produce  an  increase  tnore  than 
one-third  as  large  as  that  produced  in  model  C. 

W'e  also  e.xamined  the  effect  of  increased 
stratospheric  water  vapor  both  alone  and  in 


*Vb\.'  IvHjl  Oy  Mclimc  was  detined  as  die  ralio  of  llie  local 
O3  concentration  to  tlic  arithmetic  sum  of  the  net  chemical 
loss  rates  amt  the  bus  of  O3  through  each  level. 


Table  4.  Relative  fractions  of  integrated  odd  oxygen  destruction  rates  ascribable  to  the  Chapman,  NO^,  and 
HO^  cycles  for  various  models  over  the  indicated  domains.  For  "weighted"  integrations  the  loss 
terms  for  each  level  were  multiplied  by  the  local  concentration  and  the  local  lifetime  (lifetime 
determined  by  both  chemistry  and  transport)  ■ — ■ 


D tChapman)  D tNO'^) 1)  (IIOJ 
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1 1-45 
km 
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oombiiiation  witli  an  NO^  injoclioii  al  17  kin.  As  can 
be  seen  in  Table  2 (case  5 compared  to  case  I ),  a 
lO'l  increase  in  stratiispheric  water  vapor  tncreased  the 
column  in  model  B and  decreased  it  in  ,A.  C and 
Whett  the  increased  water  vapor  was  imposed  on 
a 17-kin  NO^  iniectioti,  the  effects  on  the  (),  column 
were  almost  additive  (case  6 compared  to  cases  2 and 

5). 

These  ex|ierimetits  have  led  us  to  the  following 
coticlusions: 

1.  A value  of  2 X lO  ' * erne's  for  Kid  leads  to 
an  reduction  due  to  a fi.xed  NO^  itijection  less  than 
half  that  predicted  during  the  ( lAP  progratn. 

2.  The  matching  of  computed  profiles  of  species 
concentration  against  observations  does  not  at  present 
provide  a definitive  criterion  for  choosing  between 
models,  even  for  models  as  different  in  their  response 
as  those  described  here.  However,  HNO^.  OH.  and 
HO,  differ  most  between  the  tnodels  and  thus  would 
provide  the  sharpest  discriminators  if  observational 
data  were  of  equal  quality  for  all  species.  The  observed 
OH  data  are  more  easily  cotisistent  with  the  lower 
value  for  Kid. 

3.  Species  concentration  profiles  must  be  known 
to  better  than  a factor  of  2 before  they  can  be  used 
to  discriminate  between  models  differing  widely  in 
sensitivity  to  NO^  injections. 


4.  The  reactions 


HO,  + 110 -11,0  + 0,. 

(KI9) 

ho"  + O3  - 116  + 20," 

(KI4) 

ho"  + HO,  - 11,0,  +0,. 

(KIb) 

110^  + no"-  no",  +" HO," 

(K2b) 

H0"+  NO,  M UNO, 

(K20) 

are  of  greatest  importance  in  determining  the  effect 
on  the  O,  column  of  an  NO^  injection  in  the  lower 
stratosphere.  None  of  these  processes  (except,  perhaps. 
K20)  has  a well-established  rale  coefficient  under 
stratospheric  conditions.  The  best  chance  of  reducing 
these  uncertainties  is  through  very  accurate  laboratory 
determinations  of  these  rale  coefficients. 

5.  Tinally,  it  is  wnrtli  noting  that  uncertainty 
resulting  from  imperfect  knowledge  of  the  rales  for 
the  Chapman  and  NO^^  desiruclion  cycles  is  small 
(less  than  a factor  of  21;  the  largest  variations  in  the 
model  results  arose  from  the  uncertainty  in  the  rales 
for  the  HO^  cycle  and  the  processes  coupling  the  NO^ 
and  110^  cycles,  it  is  also  suggested  in  Table  2 (case  .s, 
models  A and  C)  that  for  sets  of  rale  coefficients  that 
tend  to  minimize  the  response  to  NO^ , a major 
injection  of  water  vapor  in  the  lower  stratosphere 
could  significantly  reduce  the  depth  of  the  layer. 


2.2  SENSITIVITY  OF  OZONE  REDUCTION  FROM  CHLOROFLUOROMETHANES 
TO  PARAMETER  UNCERTAINTIES 


Since  the  effect  of  chloronuoroniethanes  (Cl  M's) 
on  atmospheric  o/one  was  initially  predicted  by  Molina 
and  Rowland  (1974),  a number  of  papers  have  been 
published  (for  example.  Wofsy  el  al.  IdTS,  Crut/eri 
and  Isaksen  1 975.  Turco  and  Witten  1 975 ) describing 
investigations  of  these  effects  .cade  with  various 
one-dimensional  models  of  the  stratosphere.  In  this 
study,  we  examined  the  sensitivity  of  one  such  1-D 
stratospheric  model  to  uncerlainlies  in  sonic  of  the 
imporlanl  parameters  used  in  it.  In  particular,  we 
examined  the  sensitivity  of  the  effect  of  (T'CI^  and 
CT  ,('l,  on  total  global  o/one  due  to  changes  in  key 
chemical  rates  and  diffusion  coefficients. 

The  calculations  in  this  study  were  made  with  LLL's 
one-dimensional  time-dependent  model  of  the 
troposphere  and  stratosphere.  In  this  model  92 
cheimcal  (and  photochemical)  leaclions  are  used  to 
dynamically  describe  the  slialospheric  vertical 
distributions  of  20  minor  atmospheric  species  (O('^l’). 
(),.  NO.  NO,.  N,().  UNO,.  OH.  HO,.  11,0,.  Cl.  ( I,. 

no.  no,.  (K  i().  ciNo.  cino,.  iid.  ccij‘n-,ni. 

and  CtClj)  while  three  (II.  N.  ()('l)))  are  assumed  to 


be  III  instantaneous  equilibrium.  Vertical  disiribiitions 
ot  N,.  O,.  11,0.  CHj.  and  11,  were  assumed  conslani 
thioughoul  the  calculations. 

The  chemical  rates  for  six  reactions  importanl  to 
the  chlorine  cycle  are  shown  in  Table  5.  Other 
reaction  rales  in  the  model  are  based  on  the  review 
by  Hanipson  and  Oarvin  (1975).  Solar  (lux.  absorption 
cross  sections,  and  the  calculation  of  pliotodissocuition 
rates  are  as  described  by  (iciinas  (1974). 

Sciisitivily  to  Diffusion  Coefficients 

The  vertical  transporl  in  a one-dimensional 
“globally"  averaged  inotlel  of  the  siralospheie  is 
paranieleri/.ed  through  the  so-called  eddy  diffusion 
coefficient.  There  does  not  yet  exist  a unique  approach 
for  the  derivation  of  a suitable  set  ol  eddy  diffusion 
coefficients,  and  each  nUHlelei  has  his  own  luslificalion 
for  choosing  a parliculai  diffusion  profile.  Because 
there  aie  no  uiiamhigtiotis  reasons  Idi  choosing  one 
profile  ovei  anolher.  the  sensillvity  of  the  model  to 
that  choice  nursi  be  examined.  The  various  diffusion 
coefficients  used  in  this  study  are  shown  m 1 ig.  2. 
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Table  5.  Rates  used  for  key  reactions 


Keaclion 

Rate 

Source 

Cl  + - CIO  -►  o. 

2.97  / 10  " exp(  243  f) 

/Xverage  of  Davis  e!  al. 

( 1975)  and  .Anderson  el  al. 
0975) 

CIO  + O - Cl  + (>2 

5.3  X lo  " 

^^alson  ( 1974) 

Cl  + CH^  - IlCI  + CHj 

5.4  X 10  exp(  1 (33  T) 

Davis  et  al.  ( 1975) 

on  + nci  - iijO  + 1 1 

2.0  X lo  '^  exp(-3l()  T) 

Anderson  el  al.  ( 1975) 

CIO  + NO  - Cl  -I-  NOj 

2.6  X l(f " exp(-5()  f) 

Anders*»n  el  al.  ( 1975) 

on  + nOj  - HjO  + Oj 

2 X to  " 

Low  value  quoted  hy 
lluaiipson  and  (iarvin  ( 1975) 

I'oi  each  of  tlie  lliiee  dilTusion  coclTiciciils  (('Iuiii.l' 
l‘->74.  Ciul/eii  and  Isakson  1475,  and  Wofsy  1475). 
the  following  series  of  calculations  was  carried  out; 
( 1 ) a natural  atmosphere  containing  no  chlonne.  (2)  a 
natural  atmosphere  containing  chlorine  assuming  an 
MCI  mixing  ratio  of  It)  '^  hy  volume  at  the  surface 
and  a constant  CCl^  Ilux  at  the  surface,  and 
I.M  model  2 plus  the  lluoiocarhons  Cl  (1,  and 
(T \CI,  mtiodiiced  at  the  suiface  staitmg  in  l‘>5();ind 
running  through  1475  based  on  isoild  nuorocaihoii 
production  data  (McCarthy  14^4). 

In  addition  possible  scenaiios  heyond  1475  have 
Iveii  examined  with  model  (,5a)  contmued  growth 
in  nuorocaihoii  pioduction  at  a rale  of  lO'V  per  year, 
(.'h)  continued  constant  pioduction  at  the  147.7 
pioduction  levels,  (,7c)  continued  constant  produclimi 
at  oiie-hall  of  the  147.7  production  level.  (.7d)  constant 
pioduction  at  level  until  a sharp  cutoff  in  147k 

or  m I4.S.7  (model  .7e). 

In  comparing  the  calculations  for  the  three  diffusion 
proliles.  It  was  found  that  when  compared  to  their 
resiieciive  chlorine-tree  atmospheres  the  inclusion  of 
chlonne  caused  a reduction  of  total  o/one  of  O.SCyV 
for  the  Chang  profile,  ().5o'i  for  the  Crul/en  proHle. 
and  0.(i2'-  for  the  Woisy  profile.  With  the 
chlorofhiorometlianes  (model  .7)  at  the  end  tif  1475 
there  is  a calculated  reduction  m total  o/one  (relative 
to  the  res|iective  chlorine-containing  atmospheres, 
model  2)  ol  0.77';  with  the  Chang  profile.  1.08'’V  with 
the  Ciut/en  profile,  and  1 OO'I  with  the  Wofsy  profile. 
These  results  indicate  an  uncertainly  due  to  choice  of 
eddy  diffusion  profiles  of  aiipioximalely  a factor  of 
O.k.7  to  117  Irom  the  mean  m the  o/one  reduction 
predicleil  at  the  end  of  MCTs  Inclusion  of  other 
naliiial  sources  of  stratospheric  chlorine  imiilleil  in  this 
model  would  only  allei  the  h.iseline  level  of  (),  and 
have  minimal  elfect  on  the  piedicled  pet  luihalions  due 


to  chlorolluoromelhanes,  especially  with  respect  to 
this  relative  uncertainty. 

The  differences  m predicted  effects  using  the 
different  diffusion  profiles  can  be  altrihiited  to  the 
different  rales  at  which  free  chlorine  (created  by  the 
photodissociation  of  Cl'Cl,  and  Ch,Cl,  mostly 
above  30  km)  is  transported  to  levels  below  30  km. 
The  Crul/en  and  Wofsy  profiles  allow  more  lapid 
transport  of  the  free  radical  chlorine  to  levels  where 
the  chlorine  catalytic  cycle  has  its  strongest  effect  on 
stratospheric  o/one. 

l igure  10  shows  the  reduction  in  total  o/one 
calculated  for  the  case  of  a continued  grow  th  m (T’CI, 


Year 

Fig.  10.  Sensitivity  to  the  diffusion  profile  assunuxl  is  shown 
in  these  ealeiilalions  of  the  decrease  in  total  o/one  that 
would  result  from  a continued  growth  of  lOH  per  year  in 
chlorotliioromelhane  production  beyond  1975. 


ami  ( I'ln,  priKliiction  hevitml  ol  1(J';  |hm  year 
(iiuhIoI  ,'a)  wIiilIi  coi icspomls  to  approximately  a 
sovcn-ycai  rlouhlmu  time.  As  is  expeeted.  all  preilieled 
[X'l  tiiihat''"s  increase  expoiientialls . I nr  this  scenario, 
the  rediictton  ol'  j;lo(sal  tirtal  o/onc  is  5.0  Irr  (•>.^>25''! 
by  the  eml  ol  l')‘)0  ami  1.5. P to  17. sy;  by  the  year 
2000.  It  is  sicnificani  to  note  that  there  is  only  a 
ihree-year  ilil't'erence  antonri  the  predictioirs  as  to  the 
lime  wheir  a given  reduction  in  o/one  is  reacherl. 

A constant  production  rate  at  l‘t75  levels  lor  the 
chloronuoromelhanes  (nunlel  ,5b)  was  assuinerl  in  the 
calculated  total  o/one  reductions  shown  in  I'ig.  II. 
The  shaded  area  shows  the  range  of  o/one  projections 
as  a lunclion  ot  time  when  the  I"'’.'  production  rale 
IS  used.  The  time  alter  l'i75  to  reach  a total  o/one 
depletion  ol  5'r  ranges  from  2.5  to  .5.5  years,  depending 
on  the  clu'ice  of  diffusion  prolile.  liy  the  year  2000. 
there  is  a range  of  uncertainty  ol  I ..5';'  (5. .5',?.  W'ofsy  ; 
4' ■ . Chang)  in  the  prerlicled  rerluclion  m trUal  o/one. 
Also  shown  in  t ig.  I I is  the  calculated  o/one 
depletion  resulting  fiom  a constant  production  rale  at 
one-half  the  07.5  level  (model  .5c).  At  this  lowered 
release  rale,  a meaningful  uncertainly  range  can  he 
assigned  to  the  projected  steady-state  values  as  well 
as  the  tunes  needed  to  approach  them.  The  different 
eddy  diffusion  models  lead  to  ultimate  o/one 
depletions  ol  ('.5'V  (Chang)  to  7.X"  (W'ofsy).  for  both 
scenarios  the  characteristic  limes  dilTer  very  little;  the 
time  needed  to  reach  (ine-hall  the  final  ei|uilibrium 
dcfiletion  IS  ,5fs  years  for  the  Chang  diffusion  and 
.54  yeais  for  the  Wofsy  diffusion. 


^ L_j 1 .1.1 

1980  2000  2020  2040  2060  2080 

Year 


Fig.  IT  ( aluilalerl  decrease  in  total  o/one  tliat  wrnild  result 
I'mni  a constant  prodni  tion  of  cMorotlnorometlianes  at  the 
l'#7.)  rate  and  at  one-half  the  1975  rate.  Shaded  bands  show 
range  of  \ allies  resulting  from  assumption  of  different 
diffusion  profiles. 
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The  effect  ol  slopping  pioducllon  ol  the 
chloiofluoiomcihanes  i iiown  in  I ig.  12.  I oi  this 
case  (iiunlcl  .5d)  consl.inl  pioduclion  at  07.5  levels 
was  assumed  uniil  iIk  ciul  of  l')7S  when  all  releases 
wcie  assumed  to  cease.  1 oi  each  ol  the  tested  dillusion 
piofiles.  ihe  dcsliuclion  ol  o/onc  continues  foi  some 
lime  following  the  slopping  ol  production.  I his 
happens  because  the  troposphcie,  ..cling  as  a leservou. 
continues  to  fuimsh  a source  ol  ( TCI,  andCl'iCI-.  to 
the  stratosphere.  The  miumuun  in  the  total  o/one  level 
occurs  ajiiiroximalely  at  S to  1.5  yeais  .liter 
production  is  hailed,  ilependmg  on  the  choice  ol  ihc 
diffusion  profile.  At  the  peak  reduction  in  o/one.  iherc 
was  a calculaterl  range  of  t).(r0-0.b.5'.;  addilional  o/onc 
destroyed  from  the  l'-)7b  levels,  resulting  m a I'cak 
o/one  reduction  ranging  from  1.75  to  2..527i.  /\s  seen 
in  I'lg.  12.  it  lakes  55  to  75  years  (i.e..  until  20.5.5 
or  205.5)  lo  return  to  the  already  reduced  o/onc  levels 
of  1 078  (when  production  was  halted).  I'or  the  ihiec 
dil'fusion  |irofiles  studied,  if  the  release  of 
chloronuoromelhanes  were  halted  in  1078.  o/one 
reiluction  would  continue  for  another  decade  (with  an 
uncertainty  of  i.5  years)  to  a peak  o/one  'eduction 
a factor  of  1.5  greater  than  the  1078  reduction,  and 
it  wiHild  take  55-75  years  to  return  lo  the  |078  o/one 
concentrations. 

With  production  at  107.5  levels  continued  until  the 
end  of  108.5.  total  o/one  continues  iti  fall  for 
7-10  years  with  a peak  reduction  ranging  from  2.4  to 
.5.27r'  (approximately  30'/  greater  than  ihe  108.5 
values),  and  it  would  take  .5.5-50  years  lo  return  to 


Year 

Fig.  12.  ( alculated  decrease  m Uital  o/one  that  would  result 
if  chlorofluoromclhane  production  were  slopped  in  1978 
Of  in  1983.  StiaJcJ  baiuN  show  range  of  values  ri'sulting 
Ifoni  asMimplion  of  tlifferen!  diffusion  profik*s. 


'.f 


H 


tlio  icdiiccd  I'KS.i  o/oiic  levels.  C 'onsei|ueiitl) , a 
tive-yeai  Jel;i\  in  hallmu:  produelion  nl' 

elili>ronuou>melliane  wimld  entail  an  Inerease  in  peak 
i)/one  reduelmn  by  a laetor  of  1.4  lasting  over  ati 
additional  period  of  .?0-.^0  years. 

Sensitivity  to  Keaetion  Rates 

Another  uneertainty  to  be  e.xatninerl  the 
sensitivity  to  key  reaetioti  rate  eoeffieie.it--  iinpori.  ' 
to  stratospherie  ehlorine  eheniistry.  lor  ac!.  o!  >i.\ 
reactions,  currently  used  rates  were  clianecd  t-  lest 
the  settsitivity  to  potentially  controversial  rau  choices. 
As  a reference  case  the  difUision  profile  labeled  Chang 
was  used  in  each  of  these  calculations.  Because  of  this, 
it  IS  the  ratio  of  the  chatige  in  relative  to  the 
standard  case  rather  than  the  absolute  magnitude  that 
IS  of  interest.  It  is  reasonable  to  e.xpect  that  similar 
fractional  changes  would  hold  for  the  other  diffusion 
proliles,  although  the  actual  reductions  would  be 
somewhat  different. 

For  the  reaction  Cl  + O,  — CIO  + O,.  the  currently 
used  value  in  the  model  is  ati  average  of  two  recent 
measurements  by  Davis  et  al.  (Ib75)  and  Andeisonet 
al.  (D)75).  Both  of  these  measurements  agree  quite 
well  at  stratospheric  temperatures.  However,  until 
recently  the  value  used  for  this  reaction  in  most 
models  (e.g..  in  the  model  of  Wrrfsy  et  al.  197.*^)  was 
1.85  X 10'’’.  the  value  given  in  the  Watson  (1974) 
chlorine  chemistry  review.  L'sing  this  rate  in  our 
itiodel.  we  found  that  the  reduction  in  o/one  by  the 
end  of  |9t5  increased  by  a factor  of  1.5b.  Since  the 
rates  examined  have  a ratio  of  approximately  1.7  at 
stratospheric  temperatures,  the  change  m the  o/one 
perturbation  was  nearly  proportional  to  the  increase 
in  the  rate  of  this  reaction. 

To  study  the  sensitivity  to  the  reaction  CIO  + 0 
— Cl  + (),.  we  varied  the  reaction  rate  coefficient  from 
the  value  in  Watson  ( 1974)  of  5..^  X lO"'  ' to  a value 
approximately  \l'fr  lower.  4.4  X lO  '*.  The  reason 
for  choosing  this  lower  value  is  that  it  represents  a 
recent  prelimitiary  measurement  for  this  reaction  at  the 
I'niversity  of  I’ittsburgh  (Kaufman.  private 
communication.  1975).  With  the  lower  value  the  1975 
calculated  reduction  in  o/one  was  reduced  by  a factor 


X 10'  ■ exp(- 1 1 3.5/T)  for  the  rale  coefficient  of  the 
reaction  Cl  + CH^  -•  HCl  + CHj  is  the  value  currently 
used  in  oui  model.  However,  there  is  still  uncertainty 
regaiding  this  quantity.  It  differs  from  the  earlier 
measurement  recommended  by  Watson  (1974).  5.6 
X 10'’’  exp(-l790/T).  This  rate  coefficient  hasa  large 
effect  on  the  1975  reduction  of  total  o/one.  the  old 
value  leading  to  an  o/one  reduction  about  a factor  of 
1.4  greatei  than  the  new  value  used  in  our  base  case. 
The  OH  + HCl  11,0  + Cl  reaction  rate  cociTicient 
of  2 X 10'’“  exp(-.flO.T)  measured  by  Anderson 
et  al.  (1975)  is  tlie  otic  currently  used  in  our  model. 
It  agrees  well  at  stratospheric  temperatures  with  the 
value  of  2.8  X lO"'"  exp(-410  T)  recommended  in 
Watson  (1974).  L'sing  the  Watson  rather  than  the 
.Anderson  el  al.  value  changes  the  1975  reduction  in 
o/one  by  a factor  of  0.96, 

The  reaction  OH  + HO,  — 11,0  + O,  is  extremely 
important  in  determining  the  amount  of  stratospheric- 
hydroxyl.  OH.  The  OH  distribution  likewise  is 
important  in  determining  the  conversion  rate  of  HCl 
to  atomic  chlorine  available  for  the  o/one-destroying 
catalytic  cycle.  Thougli  the  rate  coefficient  for  this 
reaction  has  not  been  directly  measured,  there  are  data 
which  support  a range  of  values  (as  recommended  in 
Hampson  and  Gai-vin  1975).  from  2 X 10'*'  to  2 X 
10'*  . Theory  and  combustion  experiments  tend  to 
support  the  lower  part  of  the  range,  although  the  low-er 
tcmjictatuic  ex^-icTiments  support  tlie  higli  end  of  tlie 
range.  If  the  range  of  values  for  this  rate  coefficient 
are  used  in  the  1-D  model,  it  is  found  that  tlie  1975 
reduction  in  total  o/one  is  reduced  bv  a factor  of  0.32 
for  the  faster  rate  of  2 X lO''**.  The  range  of 
uncertainty  for  the  reaction  Cl  + HO,  -►HCl  -i-  O, 
may  also  be  large,  but  was  not  examined  in  this  study. 
Only  the  estimated  rate  of  3 X 10'*'  cm^  s was  used 
(this  rate  coefficient  has  not  been  measured). 

Another  sensitivity  that  is  of  interest  is  the  possible 
presence  of  synergistic  effects  when  several  rate 
coefficients  are  varied.  One  test  of  this  is  to  study  the 
cumulative  uncertainty  in  the  rates  of  these  reactions 
by  choosing  the  rates  of  the  six  reactions  such  that 
they  would  maximi/e  or  minitni/.e  the  effect  on  total 
o/one.  Using  the  previously  discussed  values,  the  rates 
were  chosen  in  the  following  manner: 


•H 

The  rate  coefficient  for  the  reaction  CIO  + NO  -*  Cl 

Maximi/e 

Minimi/c 

i 

+ NO,  currently  used  in  our  model.  2.6 

Cl  -t  O, 

HI 

LO* 

i 

X 10'"*  exp(-50'T).  is  based  on  recent  measurements 

CIO  -t  0 

HI* 

10 

at  the  University  of  Pittsburgh  (F.  Kaufman,  private 

CIO  -t  NO 

LO 

HI* 

i 

1 

communication.  1975).  Until  recently,  the  rale 

Cl  + CH^ 

LO 

HI* 

coefficient  recommended  in  Watson  (1974),  1.7 

OH  + HCl 

HI* 

LO 

J 

X 10'*'.  was  used.  With  the  value  of  1 .7  X 10  * *, 

OH  -r  HO, 

LO* 

HI 

the  calculated  reduction  m total  o/one  for  the  end 
of  1975  mcreaserl  by  a factor  of  1.13, 

The  measurement  bv  Davis  el  al.  (1975)  of  5.4 


^IndiiiMcs  base  ease  rates  shtAwn  in  Table  5, 
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Ill  aiul  I ()  rel'ei  respeclivciy  to  llic  laslor  and  slower 
lalc  al  slialospheiie  lernpci alines.  .An  approximate 
laelor  ot  was  loniul  lor  the  ratio  ol  predicted  I ')75 
o/one  rediietioii  in  the  twoexlteine  cases.  This  agrees 
well  with  the  bounds  obtainable  through  direct 
nuilliphealive  superposition  ol  the  individual 

sensitivity  factors,  i.e..  ( l.56)(  1.42(1  l ie)  = 2.5  versus 
the  direct  maxiinum  laeloi  of  2..i  and 

(0.bHM0.d(>)(().5 1 ) = 0.26  versus  the  direct  ininimum 
factor  of  0.25.  C onsei|ueiilly , we  conclude  that  there 
IS  no  iinporlanl  synergistic  effect  among  these  six 
reaction  rales.  It  should  be  emphasi/ed  that  this 
inax-min  unceilamiy  range  was  obtained  by  taking 
some  exlieme  values  in  reaction  lale  coefficients  and 
does  n il  represent  the  "most  probable"  range  of 
uncerlainl>.  In  tact,  using  the  individual  sensitivity 
laclois  as  leporled  above  and  some  reasonable 
esinnales  tor  the  expected  uncertainly  in  individ  lal 
rale  coetliclenis.  one  could  derive  a "most  likely" 
uncertainty  range  that  could  be  significantly  smaller 
than  the  present  niax-min  range. 

All  ol  the  above  uncertainty  analysis  applies  to  the 
calculation  of  the  o/one  reduction  in  1675.  This 
provides  a basis  for  estimating  the  uncertainty  in  the 
steady-slate  o/one  reduction  that  is  calculated  in  a 
constant-production  scenario.  With  a constant  rate  of 
i.'jeclion  ol  (T  .Vi's  into  the  almospheie.  the  approach 
to  the  steady  stale  of  percent  o/.one  reduction  is  shown 
in  big.  II.  These  curves  are  well  described  by  an 
exponential  expression  of  the  form  .All  e‘^'(.  The 
lime-de[)endenl  term,  e‘^',  represents  the  atmo.spheric 
residence  lime.  For  the  (TM's  this  residence  lime.  1 X, 
dejx'iids  predominantly  on  the  eddy  diffusion 
coefficients  (which  determine  the  rate  of  ('FM 
injection  into  and  the  IK'I  removal  from  the 
stratosphere)  and  on  the  rate  of  ( FM  photolysis.  For 


o/one  reductions  less  than  l5'/f.  I \ depends  only 
weakly  on  the  chemic.i'  lale  constants  since  the 
amount  ol  IK'I  in  the  siialosiihere  is  nearly  unaflected 
by  changes  in  these  conslanis.  On  the  olhei  hand,  these 
rale  conslanis  sliongly  allect  the  catalytic  chain  length 
loi  o/one  desliuction.  which  aj'peais  in  the 
nuilliphealive  constant  .A.  'Hits  implies  that  a change 
III  rale  conslani  has,  to  a good  ajipioxinialion.  the 
eflecl  of  a linear  scaling  of  the  eiitiie  tune-dependent 
curve.  Thus,  the  relative  uncei  lainlies  esimialed  t'oi  the 
1975  o/one  reductions  aie  aiiplicable  to  the  steady- 
slate  limit  ol  a conslanl-iiroduction  scenaiio 

It  sluuild  be  noted  that  these  experiuients  did  not 
examine  the  lull  range  ol  uncei  tamlies  associated  with 
the  chemical  rate  inputs  to  the  model.  1 ach  of  the 
other  b(i  reactions  in  the  model  has  some  uncei l.iinty . 
and  numerous  reactions  aie  oniilled  liom  the  model 
either  because  Iheie  are  no  available  data  oi  because- 
the  reactions  are  consideied  insigiul  leant  \ 
"maximum"  and  "minimunr'  analy  sis  foi  the  complete 
reaction  set  could  lead  to  umeahslic  profiles  loi  many 
s|X’cies  in  the  ambient  almospheie.  .Such  siniulalions 
have  not  been  carried  out.  The  sensitivity  ol  the  model 
to  omitted  proces.ses  (e.g..  helerogeneous  leactions, 
ionic  reactions)  remains  at  pieseni  an  uiiitiv  stigatable 
problem,  and  pote.ilially  a .serious  souice  ol  eiioi  in 
stratospheric  models. 

■Vniong  the  reactions  that  were  studied.  Oil  ^ HO, 

11,0  -1^  O,  stands  out  as  the  most  ciiluaf  Since 
the  stratospher'c  Oil  concenttalion  is  very  sensitive  to 
the  precise  value  of  this  reaction  latc  coellicieni  and 
since  direct  laboiatory  measuremenls  ol  this  lale 
coelficieni  are  very  dilTicull.  a good  in  '.iin 
measurement  piogram  may  provide  some  very  uselul. 
although  indirect,  data  to  iiaiiow  the  lange  ol 
unceitanilies  involved. 


2.3  SOLAR  ABSORPTION  IN  A STRATOSPHERE  PERTURBED  BY  NO  INJECTION 

X 


While  much  altenlion  has  been  directed  toward 
nitrogen  dioxide  as  a chemical  reactant  capable  of 
reducing  the  o/one  layer  of  the  stialosplieie.  little 
attention  has  been  directed  toward  its  role  as  an 
absorber  ol  solai  radiation  Solar  absorption  by 
nitrogen  dioxide  had  not  been  expected  to  be 
signiltcani  conipaied  with  solai  absorption  by  o/one 
because  many  radiative  ei|uilibrinm  models  have 
generated  reasonable  stratospheric  leni)X'raUue  proliles 
without  including  nitrogen  dioxide  (e.g..  Manabe  and 
Stricklei  |9(,4).  However,  as  will  be  shown,  the 
Increase  in  solar  absoiplion  by  nilrogen  dioxide  may 
coiiiiteiact  a signiticant  fraction  ot  tlie  vlecteasc  in  solar 
absorption  by  o/one  (ollovving  sii.ilospheric  miection 
of  \()^  (oxides  ol  nitrogen). 


The  radiative  traiislei  model  used  loi  these 
calculations  assumes  a cloudless,  plane-paiallel 
almospheie  m which  there  is  moleculai  multiple 
scattering  and  gaseous  ahsor)i|ion  above  an 
isotropically  scattering  gioinul.  The  wavelength  region 
between  1X7.2  and  7.'5  nni.  in  which  molecular 
oxygen,  o/one,  and  nilrogen  dioxide  aie  the  dominant 
gaseous  absorbers,  is  divided  into  I 19  s|X'clral  intervals. 
The  atniosplierc  between  D and  55  km  is  divided  into 
4,^  layers,  which  are  I km  thick  up  to  ,55  km  and 
2.5  km  thick  from  35  to  55  km.  F.ach  almospheiic 
layer  is  divided  into  sublayers  which  have  an  opt  cal 
vlcptli,  including  scattering  and  absorption,  of  less  llian 
0.02.  There  may  be  as  many  as  500  sublayers 
ik'pending  upon  the  total  optical  de|'lh  ot  the 


alimisplioio  llio  (i;u:s!.-Sculi.-I  lleialivf  sclicino  is  used 
111  solve  the  ladialive  liaiisl'er  ei|uatioii  lo  oblalii  the 
radiative  inlensity  at  mcienieiits  ot  h°  m the  local 
/eiiith  anple  at  each  level,  ( oiiipoiients  of  the  radiative 
intensity  aie  then  inteitrated  ovei  the  upper  and  lowei 
lieniispheres  to  obtain  the  diH'use  tluxes. 

I he  absorption  cross  sections  were  derived  Iroin  a 
variety  ol  sources  described  in  (lelinas  (l'*74).  The 
vertical  pioliles  ol  teinperature,  pressure,  and  oxygen 
concentration  correspond  lo  the  IhS.  Standard 

■Alinospheie  (I‘>h2).  and  the  concentration  protiles  lor 
(),  and  NO,  were  derived  I'roin  Iransport-kinetics 
calculations. 

I he  two  cases  considered  are  NO^  iniections  (as 
NO,  I at  I"  or  20  kin  at  the  rale  of  2000 
molecules  enr  -s  lUiil'orinl'  distributed  over  a 

l-km-lhick  laser.  This  injection  rale  (2..x 
I * 

X 10  ■ g year  I corresponds  'o  that  of  a Heel  ol 
several  thousand  suixtrsontc  transports  ((Irohecker 
el  al.  The  stratospheric  column  densities  ol  Oj 

and  NO,  lor  each  ot  the  three  cases  are  shown  in 
Fable  h.  The  tintx’i tiirhed  O,  profile  corresponds  to 


Table  6.  Stratospheric  column  densities  (above  13 


km) 

for  0.,  and  NO, 

(molecules/cm  ) 

(a>c 

'>3 

\() 

2 

1 nperlurhcd 

7.92  X lo'** 

4.06  X 

lo'-' 

17-km  injection 

7.50  X I0‘’^ 

5.43  X 

to'-^ 

20km  injection 

7.03  \ 10'** 

6.80  X 

lo'-*’ 

0..'().s  ainrem  ol  ()^.  The  I '-km  injection  corresponds 
(by  ntodel  calculation  I lo  a rciluclion  in 

sliatospheric  ().,  and  a .'3.7'.'  increase  in  slralospheric 
NO,,  whereas  ihe  20-kin  injection  corresponds  lo  an 
I 1. 2'.'  reduclii.)n  in  O,  and  a b7..S'7  ineiease  in  NO,. 
In  each  case,  the  fractional  change  in  ihe  NO,  column 
density  is  approximately  six  times  the  fractional 
change  in  the  O,  column  density. 

For  each  of  the  three  sets  of  concenlralion  profiles. 
O^  and  NO,  solar  absorption  rates  were  computed  for 
values  of  the  solar  /enilb  angle  (I  langing  from  0°  to 
7.x  aiul  values  of  the  surlace  albedo  a ranging  fioin 
0 to  I.  For  the  un|H'rlurbeil  case,  ihe  siralosplrei ic 
O,  absorption  langed  Irom  2.'t.s'.  ( d ~ 0 . a - Ol 
lo  ^.23'-  (0  7X  . a._  ■ I)  ot  the  incoming  solai 
ladialion.  For  the  unpeiluibeil  case,  the  slialosplieiic 
NO,  absorption  ranged  trom  0.072''  U>  - 0 . a_^  ' 01 


lo  0.365';  (d  = 7X  . a^  = I I ol  the  solai  ladialion 
incident  al  the  top  ol  the  atmos(iheie.  Fills  result 
verities  that  ihe  solai  absoipllon  bv  NO,  is  noimallv 
only  a small  fraclion  (0.024  to  0.050|  ol  the  solai 
absorption  by  Oj. 

The  ratio  ol  the  increase  in  siralospheiic  solai 
absorption  by  NO,  (with  respect  lo  the  unpei imbed 
case)  lo  the  decrease  in  absoipllon  by  O,  is  given  in 

Table  7.  The  ratio  of  the  increase  in  NO2  absorption 
to  the  decrease  in  absorption  for  a 
surface  albedo  of  0.26  '•  ' 

Solar  /enilh  angle 

Case  0 30  60  7X 

1 7 km  injeelion  0.38  0.40  0.45  0.50 

20-km  injection  0.35  0.37  0.41  0.47 


Table  ~ toi  selected  solar  /enilh  angles  and  a surface 
albedo  of  0.25.  This  ratio  varies  significaiitlv  with  solai 
/eiiith  angle  but  only  slighllv  with  surface  albedo.  The 
ranges  ol  this  ratio  lor  the  I”-  and  20-kni  injeclions 
iiulicale  that  Ihe  increase  in  NO,  absorption 
significanllv  compeiisales  (30-50'.')  for  the  decrease  in 
().,  absorption. 

The  change  111  sliatospheric  composition  increases 
the  solar  radiation  iransmilted  lo  Ihe  lioposphere. 
resiilling  in  increased  tropospheric  heating.  More 
radiation  is  scattered  upward  from  the  Iroposi'liere  and 
a laigei  fraclion  of  this  ladialion  escapes  to  space 
iveause  of  the  increased  transmission  of  the 
siralosphere.  'Fins  leads  to  an  increase  in  Ihe  planetary 
albedo  and  a decrease  in  Ihe  net  solar  heating  of  the 
alniospheie  earth  system.  .\n  analysis  ol  the  perluibed 
solar  radiation  budget  lor  the  20-km  inieclioii  case  is 
shown  111  F ig.  13  foi  a solar  /enilh  angle  of  60  . The 
increase  in  sliatospheric  absorption  by  NO,,  tbe 
increase  in  (roposplteric  healing,  and  the  increase  in 
radiation  lost  lo  space  combine  lo  equal  the  decrease 
in  stiatospheric  absoipiion  by  (),.  The  lelalive  niag- 
niliides  of  these  various  componenls  van  w ith  surface 
albedo.  F oi  small  values  ol  surlace  albedo,  most  of  the 
net  decrease  in  slralospheric  liealins;  goes  into  in- 
creased iropospheric  healing.  I'oi  large  values  of 
surlace  albedo,  iiiosl  of  the  net  decrease  in 
siralospheri,  heating  is  lost  lo  the  syslein  as  radialion 
scattered  lo  space.  Ihe  increase  in  NO,  absoipllon  is 
aiiirioximaiely  a coiislant  fraclion  of  the  decrease  in 
(),  absorption  foi  all  values  ol  suilace  albedo,  as  was 
piev  loiisly  men  honed. 
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Fig.  13.  C'alculatod  decrease  in  sirulospheric  absorption  by 
(upper  curve)  and  the  three  components  that  make  up 
(his  c|uanti(>  as  a tunction  o(  surface  albedo.  Ihis  evample 
IS  for  the  2(kkm  mieclUM)  case  and  a solar  zenith  angle  of 
bO  . 


■Mlhoiii;li  [XTUirhalums  lo  the  sliatosplierie 
eompi'Sition  due  to  \()^  mieellon  letul  to  inereased 
solar  healiML’  ol  the  iroposf'iiere.  it  does  not  iioeessarih 
tvillovs  that  these  pertiirhatioiis  eaiise  an  increase  in 
surtace  tent(X'iaUirc.  I’evtvirhalioirs  to  the  loixtsvase 
(IK)  ladiation  hudeet  must  also  he  considered  wiieii 
computing  clianees  in  surl'ace  leiii|X’ralure.  Because  ol 
the  lediiction  in  solar  heatine  ol  llie  stralosphere.  the 
slialospheric  temperature  smII  decrease.  I his  reduction 
in  emittme  temperatuie  comhined  with  the  rcdticUon 
111  o/one  coiicenli.it:  111  results  in  a decrease  in  ihe 
iliuMiuaid  lonevtave  radiative  I'lu.c  eiiterinu  the 


tiopospherc.  lliis  cooliny  elTect  lends  to  counteract 
the  increased  solar  heating  ol'  ihe  tiopospherc. 
Kauianalhan  cl  al.  (I‘I75)  found  llial  the  longwave 
el  feel  dominales  ihe  solar  efieci  for  glohal  mean 
condilions.  I'sing  a radialive-conveclive  model  and  a 
ix'rtuihation  similai  to  the  dO-km  mieclioii  case,  lhe> 
computed  a change  in  sin  lace  temperatuie  ol  0.1  K 

.Al  hiologicallv  nn|H)ilanl  wateicngllis 
(2SO-3i^O  mn).  solar  ahsoi|nion  hv  NO,  has  only 
a small  effect  on  the  amouni  ol  solar  ladialion  leaching 
the  earth's  surlace.  In  this  six-cttal  region  the  optical 
thickness  ol  {).,  is  seveial  oiders  of  magnilude  giealer 
than  the  optical  thickness  ol  NO,;  Ihus  ihe  llu.s  of 
l'\  -B  radiation  incideiii  al  llie  eailh's  surlace  is  much 
more  sensitive  lo  changes  in  O,  column  deiisily. 
Neglecting  solar  ahsorplion  hy  NO,  leads  to  a shglit 
overestimate  of  the  hiological  effects  from  slialo- 
spheric  injeclioii  of  NO  . I'or  tlie  17-  and  ,!!0-kni 
inieclion  cases  considered,  iieglecliiig  solar  absorplion 
by  NO,  caused  Ihe  change  in  radiation  at  Ihe 

earth's  surface  to  be  overesliiiiated  by  Z'  . 

The  fact  that  a sigiiificaiil  fiaction  of  the  decrease 
in  net  stralosphei ic  healing  is  lost  to  ilie  sysleni  as 
railialion  scattered  to  s|iace  has  uiiportaut  implications 
foi  radiative  er|uilibiium  models  and  clinialic  moilels 
depeiideiit  on  radiative  ciiergv  balance  calciilalions. 
The  effect  on  planelarv  albetlo  of  changes  in 
atmospheric  composition  must  be  included  to 
accuratelv'  comptile  the  peiiuibed  .solai  radiaiion 
budget.  Noi  including  the  albedo  effeci  nia>  even  lead 
lo  predicted  changes  in  surface  lempei alure  r'f  the 
wrong  sign  (Ramanaihan  et  al.  IslTsl.  The  effect  of 
solai  absorplion  b\  NO,  should  also  be  included  in 
lliesc  models  loi  case  studies  involving  signilicani 
changes  in  NO,  column  densilv  . Othei  species,  such 
as  water  v.ipoi  and  aerosols.  shouKl  also  be  inchuled 
ill  more  coiii|ileie  anaivses  of  SSI  effects,  but  the> 
woiikl  not  affect  the  c.|ualualive  results  piesenied  here. 


2.4  EFFECT  OF  NO  PHOTOLYSIS  ON  NO  MIXING  RATIOS' 


It  has  been  suggested  that  the  oiiK  significant  sink 
lor  odd  nilrogen  in  the  stratosphere  is  ihrougli 
iranspoii  ol  UNO,  lo  the  trofuisphere.  As  a result, 
it  NO  were  to  be  miecled  m the  middle  or  lower 
slr.ilosphete.  the  change  in  odd  iiiltogeii  INO^  “ NO 
+ NO,  ♦ UNO,  + N|  mi\,ing  ratio  above  the  height 
ol  mieclion  would  he  ne.irlv  corislanl  (National 
As.idemv  ol  Sciences  l'i'5b.  pp  II4  II‘0. 

I’hotoKsis  ol  milk  oxide  and  the  siibsetiiiin I 
reactions  ol  imrogen  aloms  with  nitric  oxule  act  as 


a sink  fill  odd  niliogen  in  the  uppei  stralosphere.  This 
papei  lepoiis  an  investigation  of  ihe  iniporlaiice  of  Ihis 
net  \0^  sink  by  means  of  a one-diiiiensional 
Iransport-kmelics  model  (('hang  el  al.  I')74)  hi  Ihis 
model  the  leaciions 


hr 


NO  • N + O. 

N + NO  -•  N,  r O. 

N V No,  - N,0  + O 
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piovkk'  a dicniical  sink  lor  inld  iiilmuL’n  in  tin.'  n|'ix-i 
siialosplKMC.  I'sinp  ilio  C'lianp  (l‘)74)  eddy  dilTnsion 
proiile.  a paianictcii/aluin  ol  NO  pliololysis  as  shown 
b\  Cieslik  and  Nieolel  (!'»?.').  and  ilie  elieinieal 
leaetions  used  h\  Dnewerei  al  ( l'>70).  «e  earned  onl 
tire  lolhnv  u;a  es.penn\eius; 

1.  We  eoinpnled  an  anihienl  N()^  iniMnu  lalio 
and  a peiinibed  NO^  inivinu  laiio  due  lo  an  inieetion 
ot  NO  at  a rate  ol  2000  nioleenles/enr^-s  in 
a l-kni-ilnek  layer  eenieied  at  l~  kin  iisina  onr 
Niandaid  eheinislrs  (wiili  the  Kite  eonstant  I'oi  Oil 
r no,  - II, O + O,  ecpial  to  2 X lO  " ein^  We 
then  repealed  the  ealeulation  with  tlie  rales  ol  the  NO 
pholoKsis  sink  set  eipial  to  /eio. 

2.  both  ol  the  above  sets  ol  ealenlaltons  were 
reiieated  Ibi  the  Ihinten  etld>  dilliision  proiile 
(Cliniatie  lni|iael  ('oinniiltee  lO's,  p,  42)  and  lor  the 
('nit/en  eddy  dilTtision  profile  l('iiii/en  and  Isaksen 
l»?5i. 

lire  NO^  inisinr;  latios  resultiin:  Iron)  tlie  abtwe 
ealeiilations  are  yiren  in  I'lt!.  I4.  A eoini'arison  of  the 
lesnlts  with  and  without  the  NO  photolysis  sink 
iniiihes  that  this  sink  lor  NO^  is  lesponsiblc  lor  holh 
deeieasinu  [NO^  | al  all  altitudes  and  lot  a decrease 
in  NO^,  niiMiip  ratio  with  altitude  above  about  oO  km. 

It  IS  worth  noline,  liowevei.  that  for  the  peiltirberl 
eases  the  ehanees  in  iniMiie  ratio  are  not  eonstant  with 
alliliule  even  when  the  phololvsis  sink  I the  onlv 
sheinieal  sink  Tor  NO^  in  the  model)  is  removed.  This 
IS  a lesiill  ol  a ehanee  in  the  net  rate  of  the  "naliirar' 
proihietion  of  \0^  from  the  reaetnm 

N,0  + Ol  '|)l  • 2NO. 

This  somee  l'■lnl  ehanees  heeanse  the  o/oiie 
eoneenliaiion  changes.  The  o/one  redtielion  leads  li> 
a decrease  m Ol’l))  eoncentraiions  al  high  allilndes 
where  the  pioportional  change  nr  sotai  lliix  is  less  ihan 
Ihe  ehanee  in  o/,one  eoneeniiation.  and  an  ineiease  m 
0(  'l)|  at  low  altiliules  where  the  ehanee  in  solar  iTiis 
density  more  than  eonipensales  for  Ihe  ehanee  in 
o/one.  The  ehanees  in  N,()  Ibeeanse  ol  altered  solar 
Ihi.x)  aie  relatively  niinoi.  The  net  change  in  the 
production  ol  NO^  liom  0('l))  + N,0  m the  model 
will)  NO  photolysis  suppressed  is  given  m Tig.  Is  foi 
the  three  diffusion  profiles. 

Ihe  net  changes  m Ihe  o/one  eoliimn  foi  Ihe  above 
ealciilatioiis  are  given  m Table  X.  In  e\amining  the 
o/one  rerinciions  il  should  be  noted  that  the  rate 
coellicienis  used  in  the  current  c.ilculalions  aie  the 
recommended  values  liom  llampson  and  (iaiv'in 
(|v)75),  with  . X lo  ''  ciii'^ 's  for  Oil  r HO,  Thus, 
there  aie  numeo'us  minoi  changes  and  one  niaior 
shanev  in  lire  Kile  ciK-llicienis  used  here  comparevl  to 
those  userl  to  eener.ile  ihe  ( I \l’  findings  (< iiobecker 


Fig.  14.  IMiiis  ot‘  ,aluiiatoil  amlrieirl  NO  niistng  rjito 
Kiinliienl)  and  changes  in  M)  niising  ratio  (A.V)  I alter 
a 17-kin  ,NO^  intection  as  descritred  in  the  test.  .Solid  curves 
are  calcidalions  willi  Itie  normal  ntiHlel.  dashed  curves  are 
calculations  with  the  NO  |)hololysis  sink  suppres.sed. 


el  al.  l'>~4).  Allhough  omission  of  Ihe  N ati'in  sink 

foi  NO^  makes  a veiy  substantial  thffeience  in  the 

NO  inKiiie  ralui  above  .7(1  km.  it  causes  only  a minoi 
> 

incre.ise  m the  percentage  o/one  pei lurbalioti  lor  NO^ 
inieclions  I Ins  results  because  both  ilie  ambient  and 
iviluibeil  NO^  mismg  ratios  .ne  reduced  by  ihe  sink. 
The  o/one  peilurbalion  is  mote  slionglv  alieeieii  fot 
the  llunten  profile,  which  was  used  in  ihe  foimula 
developed  by  NAS  (Naiional  Academy  ol  Sciences 
HI7,sh.  p,  d.’l.  than  loi  the  olhei  di'itision  pioliles 
ei'iisideied. 
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It  IS  iiidiiiJcil.  A iiukIcI  llul  a)iii|Hiicd  a U.51'' 
rodiK'lluii  111  o/oiK-  diic  tvi  diloiutluuios'ai lions  111  1'I7(' 
usiiis;  oiii  till!  clioiiiisii)  aniipiiknl  only  a O.dU' . o/onc 
icdiiclion  111  l't"’(i  uIkmi  no  pliololysis  was  iicglcckHl. 
1 Ills  I’llik'l  oocilis  bsVaiisi.'  llio  kiaclioii 

NO  1 ( lO  • NO,  + Cl 

kklik'cs  llio  crtk'ioiicy  of  both  the  chlorine  and  NO^ 


o/oiie  destiiielioii  e>eles.  and  ihe  elleet  ol  this 
leaction  is  es|veiall\  de|ieiideiit  on  the  NO^  iiiixiiij! 
lalio  above  .'.s  kin. 

I Inis,  allhoiieh  models  ol  the  ellecis  of  NO 
mieclioiis  on  the  slialospheie  aie  not  \ei\  seiiMlne 
to  ihe  NO  photolysis  sink  loi  NO  . inclusion  ol  this 
sink  Is  i|iilte  iiiipoilaiil  loi  iiioilels  of  the  efiects  ol 
chloiollnoiocai bolls  in  the  stiatospheie.  Indeed,  ihe 
concentiatioii  ol  NO  above  .'li  kin  is  an  iinportanl  and 
)Hioily  valulated  i|iiaiitity  in  cm  lent  models. 


2.5  OTHER  INVESTIGATIONS  USING  THE  1-D  TRANSPORT-KINETICS  MODEL 


Traiisieiil  Effects  in  the  Transport-Kiiielies 
Model 

When  the  tiiiie-depeiuletit  response  of  the  l.Ll. 
one-diniensional  model  to  a cotrslant  NO^  injection  is 
computed,  there  is  an  initial  period  when  the  model 
computes  an  increase  or  very  slowly  growing  tlecrease 
in  total  o/one.  .An  initial  increase  in  total  o/onc  is 
obtained  when  the  rate  of  the  reaction  of  (.ill  with 
HO,  is  taken  to  be  near  2 X 10''  ' enr'  s (model  A. 
desciibed  eailier  in  .Sec.  2.1  ) and  the  in|eclioii  alliliide 
IS  I km  Ol  less. 

This  behavior  arises  because,  m the  region  below 
about  20  km,  a local  increase  in  NO^  results  iti  an 
increase  in  the  local  conceiilvalion  of  o/vuie  in  lliat 
model,  and  if  the  iniectioii  altitude  is  iiuich  below 
20  km.  o/oiie  is  increased,  until  injected  NO^  is 
traiisjioi ted  up  into  the  legioii  above  20  km  where  it 
causes  o/one  destiiiciion. 

Ihe  increase  in  o/one  at  lowei  alliliules  results  lioni 
a combination  ol  thiee  laclors.  Ciiveii  the  rate 
coeflicieiil  Used,  the  iiiosi  inipoiianl  ol  these  (in  terms 
ol  the  total  o/oiie  coluniiil  is  the  mlerleieiice  m IIO^ 
cycle  viesliuclion  ol  o/one  ..aiiseil  by  Ihe 

inlunalely  related  ineihaiie  smog  leaclion  has  a snnilai 
elTecl.  .Old  is  ihe  inaioi  somce  ol  the  Iropospheiic 
o/one  mciease  Ihe  elleci  ol  high  altitude  lediiellons 
111  o/oue  on  liie  ladiatioii  reachiii!',  Ihe  lovvei 
stralospheie  results  in  ,ui  mciease  in  lovvei 

slialospiieiic  o/one  at  lalei  limes,  but  h.is  little  elfeci 
on  Ihe  eaily  liaiisient  hehavioi 

Ihe  I I I irealmeiit  ol  methane  ovidaliou  uses  the 
lollowmg  leaclions 

( 11^  • OH  - 11,0  + ( II,. 

( 11^  ^ 0('l)|  •'‘oil  t ( ii, 

(),.  NO 

( 11-  ■ • 2HO,  r (.  O (assiiiiied  to  •,ccur 

lapidly  1. 

Hits  slioilciii  mechauism  yields  loughly  the  coned 
imomil  ol  HO,  |iioduetion  while  mideiesimialing  O, 


production  by  about  oiie-lliird.  It  avoids  an  mciease 
in  the  number  of  species,  and  piodiices  sigiiilicani 
errors  only  in  the  troposphere,  which  is  controlled  by 
boundary  conditions  in  any  case. 

The  transient  bcliavioi  calculated  using  the  ( hang 
(HI74)  eddy  dilfusion  profile  (f  ig.  2l  foi  a 1 7-kni 
inieclion  of  NO^  ai  a rale  of  2000  molecules  enr  -s 
in  a I -kill-thick  layer  is  illusiraled  m Tig.  lb,  and  the 
iransieiil  behavioi  for  a lO-kiii  NO^  ui|eclion  ofevpial 
magniliide  is  illusiraled  in  Tig.  1 7.  Tlie  references  lo 
hianch  1 oi  blanch  2 concern  the  assumed  |iioducts 
of  NO,  photolysis.  In  calculating  the  curves  labeled 
branch  1 . NO,  photolysis  was  assumed  lo  y held  NO 
•I-  (),.  aiul  in  calciilaiing  those  labeled  branch  2.  NO, 


CM 


CO 


o Time  - years 

Fig.  16.  ( alciilaleil  cliaiige  in  Ilie  lolal  O,  coliinin  vs  lime 
as  a resvill  ot  an  mievlion  (as  NU,1  ol  200(1 

moleuiles  enV  -s  in  a I kni-llmk  layer  eenlered  at  17  l,m. 
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AO^  (total  column)  - 10  ^ molecules/cm‘^ 


Time  - years 


Fig.  17.  ( alculatcd  cltjiigc  in  llic  (o(al  ct»!umn  vs  time 
as  a n-Mill  of  an  NO^  injcttion  (as  NO,)  of  2000 
inolsiiik-s inr -s  in  a I Kni  lhii  k layiT  scnfcri'il  a(  10  km. 


^ ()  well'  ihc  assiiiiK'il  |)iodiii'ls.  ( miia'  IX  pifseniN 
dll'  alliliuli'-ik'i'i'iuleiil  I'liaiific  iii  o/oiii'  i-oiUi'iiOalioii 
aliei  U.2.  U.S.  2.  5.  ami  11)  veais  loi  llie  lO-kni  NO 

s 

iiili'ilioii  assiimiiii'  hiaiii'li  1 

I III'  ii’Mills  ili'Si  iihcil  Ill'll’  ail'  111  qiialilalivi' 
Mjaiilii'aiiii'.  I lii'V  all'  si’iisilui'  lo  llic  ilion'i’ ol  pooiK 
known  lali'  lonsianis.  .iiul  aie  i’\pi’ili'il  lo  Ik-  ipiik' 
si’iisilivi'  lo  till’  I'llils  ilillusion  piot'ili'.  Till'  tneali'sl 
uiii'i’ilainu  IS  assoi'ialfd  willi  low  alliuid.'  iniCiUons 
ol  NO^.  I inis,  ilii’  i|nalilalivi'  a'.snlis  slioiild  be  iiealed 
lanlioiisl) . 

New  Eddy  Diffusion  CiK'fficieiit  Profile 

new  vertical  eddy  diltiisioii  eoetlieieul  piofile  has 
been  deiived  loi  ibc  M ) lianspoi t-kinelies  model.  Tins 
pi'ollle  was  designed  lo  improve  Ibe  fil  wilb  inetliane 
and  nuawis  oxide  ineasuicinenls  in  the  ripper 
sliaiospliea’.  Il  was  neeossilaled  by  llie  a'ljniri’ini'ni 
lor  im[)ioved  lepresenialion  of  liaiisporl  in 
llie  ripper  siralospliere  where  the  distiibrition  ol 
ehlorol Irioronu’lhanes  beeomcs  impoilanl.  dlie  new 
I'drK  dillrision  eoeflieieirl  profile  (shown  in  Tig.  1')) 
provides  loi  inerea.scd  lianspoi  l in  ihe  ripper 
slialosplieie. 


9 3 

AO,  - 10  molecules/cm 

o 


Fill  18.  (.iluilalid  ihaiijie  m o/one  r imir'iilralioii  vs  aliiiiiile  after  0.2.  O.S.  2.  .S.  and  10  years  (<ir  (Ire  Kkkm  Ml  inieelion 
»t  I 1)1  1 ' 


A1 ti tude 


fcl'fect  of  CIONO,  in  the  Stratospheric  Model 

I he  el  led  ol  addiiii:  CIONOt  io  tlie  l -l) 

iranspoi l•klnclies  model  is  shown  in  Table  9.  The 


50  f--  Chang  ( 1 974), 
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Chang  (1976) 


Vertical  eddy  diffusion 
2 

coefficients  — cm  /s 

Fig.  19.  New  vertieal  eddy  diffusion  eoeffieicnl  profile 
(Chang  1976)  eonipared  with  the  previously  used  profile 
(Chang  1974). 


adiial  leaetioii  rate  for  lorination  of  ( l()\0,  is 
unknown  at  present,  but  preliminary  measurements 
(Kaufman  1976.  Birks  I97(i.  IX'More  I97(r.  ail  private 
eommumeation ) hav  ■ mdiealed  the  rale  coefficient  is 
between  O.O.S  and  'J.l  times  that  lot  loitnation  of 
UNO,.  Case  0 refers  to  the  model  vrithout  CIONO,. 
In  case  1 the  ( lONO,  productioti  rale  was  assumed 
to  be  O.i  tunes  that  of  IINO3.  and  in  case  2 the 
CIONO,  prodtiction  rate  was  assumed  to  be  0.05  times 
tluit  of  UNO,.  A reaction  rate  of  2 X 10"'^  cm^  s 
(Oasis  |97(i,  private  communication ) was  used  lot  the 
{ lONO,  + 0(3p)  in  cases  I and  2.  The  photolysis  rate 
calculated  111  the  model  for  CIONO,  includes  the  rc- 
ceiitly  revised  absoiption  cross  sections  from  .Spencer. 
Molitia.  and  Kowland  (I97(s).  and  the  new  Chang 
ditiusion  piofile  I Tig.  19)  was  used  m all  calculations. 

As  seen  m Table  9.  the  effect  of  CIONO,  is  to 
reduce  the  effectiveness  of  the  CIO^  released  by  the 
chloninuoromethanes  by  approximately  a factor  of  2. 
With  constant  prodiiclioii  of  the  chloroHuoromethaiies 
at  |97.s  levels,  the  o/one  reduction  by  the  year  217b 
changes  from  without  CIONO,  to  7.47'  ' for 

case  2 which  has  CIONO,  produced  at  0.05  times  the 
rate  of  IINO3.  In  the  model,  the  effect  of  CIONO, 


T-.fele  9.  Percent  change  in  total  ozone  from  ambient  atmosphere  (containing  CCI^  and  CH^CI)  due  to 
chlorofluoromethanes  — 


t ase  0 (vvithoul  CIONn^j 

Case  I (with  nONO,  = O.I  UNO,) 

Case  2 (with  CIONO^  = O.O.S  IINO3) 

Case  2 (with  CIONO  3 and  rate  change): 

Oil  + HOj  ■ lIjO  + (>2 
(2  X l(f  ’ 2 X l(f"S 

(I  + O,  • CIO  + Oj 

(2.97  X l(f " 1.8a  x lo  ' ’) 

no  + o ■ n + Oj 

(3.38  X l(f " e^'*’^;  5.3  X l(f ") 

no  + NO  ■ n + N(>2 

(1.13  X l(f ' ' e^"”  1.7  X lo  ' *1 

( I + ni^  iin  * ni, 

(5.4  y Kf'-’  e "^^  ' : 5.6  x l(f"  e '^'^"  '1 

OH  + iin  ii^o  + n 

(2.0  X l(f'^  e^'*’  2.8  X l(f'^  e”^"’  ') 


1974 

1976 

2i7f»  (constant  production 
al  1973  levels) 

0.61 

0.86 

13.68 

0.23 

0.30 

4.87 

0.38 

0.5 1 

7.47 

0.02 

0.03 

0.43 

0.59 

0.38 

0.53 

0.40 

0.55 

0.41 

0.56 

0.33 

0.45 

Concentration  - niolecuies/cm 


Ill  CIOSO,  tor  holli  c;iM's  1 ;iik1  2 arc  less  than  llic 
ii|'l'ci  Imiii  ot  lo'*  inolccules  ciir'  civeii  by  Mincray 
( I‘)7m. 

Also  slioun  ill  Tabic  '■>  is  ibc  cUccl  of  varyinc  ibc 
same  si\  icactloiis  discusscil  above  when  CIONO,  is 
melnilei.1  m the  model.  Tiie  rales  shown  m I'aienlheses 
indicate  the  rates  noimalK  used  in  the  model  Ion  the 
left)  and  the  rales  used  to  lest  the  iiHidel's  sensitiMiy 
(on  the  right).  Whereas  ehanging  the  lale  eoelTieienl 
loi  on  + no,  liom  : X lo  " to  : X lo  "'  em^-s 
had  lerlneed  the  elTeeliveness  of  (10^  In  a laeloi  ol 
,'.l  wiihoni  CIONO,.  die  effeel  with  CIONO,  is  to 
reduce  the  change  in  total  o/one  hy  a tactoi  ot  17 
(Ironi  0.-^1'.  lo  O.O.V  ) foi  the  year  l‘>7(i.  This 
increased  sensitiMlv  is  atlribriled  lo  CIONO,  l>ing  up 
NO^  in  an  atmos[iliere  wheie  O^  is  predominaiitlv 
controlled  In  NO  when  Oil  + HO,  is  2 X lU"'*' 
enr  s (this  rale  will  lediice  the  Oil  available  lo  react 
with  NO,  lo  produce  UNO,,  lesiihmg  in  increased 
NO  and  NO,  levels  lelalive  to  when  die  late  is 
2 X 10'").' 

Diurnal  calculations  have  also  been  carried  out  foi 
the  model  containing  CIONO,  Pie  dim  mil  variations 
r'f  CIO  and  CIONO,  afler  the  I2lh  day  (staitmg  from 
tlie  constant-sim  nurdel  results)  aie  shown  in  l ig.  20. 
Tlie  diurnal  calculation  shows  less  dinrnally  averaged 
CIONO,  as  compared  with  values  liom  the 
consianl-siin  calciilalion  (see  crosses  in  f ig.  20).  This 
is  due  lo  conversion  of  some  ol  the  CIONO,  lo  PCI 
througli  the  diurnal  variation  ol  NO^  and  CTO^  species. 

Effects  of  O,  ^ 0('D)  + O,  Quantum  Yields  Near 
Threshold 

In  comparing  our  model  results  wiih  similai  lestills 
front  the  Michigan  and  NC.AK  gioups.  we  lotiiid  that 


die  III  model  v lelded  aboiil  hall  the  total 
slialospheric  NO^  contained  in  the  olhei  models  Ihe 
reason  was  that  the  other  modeline  groups  weie 
assuming  a unit  nuannim  yield  toi  0(^  0)  production 
Irom  o/one  photolysis  al  wavelengths  shoilei  than 
.Ot)  iini  and  /eio  al  longer  wavelengths  (as 
lecommended  by  llamfison  and  Oarvin  I'I75  ).  whereas 
we  were  using  the  lornmiation 

Oi'l),  = I I l"i  X < .t05  mn. 

I + e'*  ' *’ 

I 

"I'l'l  = lai  .tOS  < A < .t.tO  nm. 

I , ^.lA  .ULt)  4 S 

which  was  recommended  by  Johnston  (P)7.v).  I se  ol 
Johnston's  expiession  reduces  the  lower  stratospheric 
0('l))  production  by  about  a taclor  ol  2 relative  lo 
die  .'lO-iini  step  lunction  and  accounts  for  most  of 
the  differences  in  NO^  between  the  l.LL  model  and 
the  othei  models.  Because  an  abrupt  threshold  is 
theorelically  unlikely  and  seems  incompalible  with 
recent  results,  we  believe  the  form  of  e.vpression  used 
in  our  model  to  be  su|ieiioi  to  the  .MO-nm  step 
function.  However,  we  recogni/e  that  the  temperature 
dependence  and  precise  value  rd  0('D)  in  die 
.JOO-to-.i  I .‘'-mil  range  is  i|tiiie  uncertain  Tims  the 
expression  used  is  signil'icantly  unceilain.  The 
sensitivity  of  stratospheric  NO^  to  an  incieniental  NO^ 
miection  is  fairly  sensitive  to  the  ambient  NO^ 
concentration  (in  an  atmosphere  with  doubled  NO, 
die  response  of  total  o/one  m model  A lo  a l7-kni 
NO  miection  changed  from  a base  case  1 
decrease  to  a I ..'.s'  - decrease).  .Also,  depending  on  the 
iiupoiiance  ol  chlorine  nitrate,  the  response  rd  o/oiie 
lo  a source  of  slialospheric  chlorine  may  be  rpiite 
sensitive  lo  the  anibienl  NO^  conceniration. 


2.6  COMPARISON  OF  MODEL  RESULTS  WITH  OBSERVATIONS 


Catalytic  Destruction  of  O/one 

We  have  been  collecting  observational  data  believed 
lo  have  a beaiing  on  the  validity  of  the  catalytic 
destiuclion  ot  sliatospheiic  (),.  Pie  data  i|uahlalivelv 
conlirm  the  current  theory  foi  shorl-peiiod  responses 
ol  O,  al  levels  above  .'()  km  lo  ('olar  cap  absorption 
events  (I’t  \'s)  and  changes  in  solar  fins.  Above  .'()  km 
the  strati  -I'heie  is  m photochemical  ei|uihhinim.  and 
o/one  should  be  contiolled  by  ['holocheniislry  in  this 
region  The  data  reveal  that  dynamics  may  be  the 
dominant  controlling  ['rocess  allecling  long-leim 
trends  in  o/one  below  .JO  km. 

Specdic  examples  ol  dilleiences  between  model 


results  and  observational  data  are  cilcd  below.  I hesc 
differences  reveal  aieas  wheie  we  have  limited 
understanding  about  the  piocesses  dial  may  be 
alfecliiig  o/one.  I uilher  reseaich  is  needed  lo 
satistaclorily  explain  these  dd'Icrenccs. 

1.  O,  and  nucleai  tests.  .Since  l‘i('2.  loial  (),  has 
incieased  to  mean  levels  liighei  than  obseivcd  beloie 
atmospheric  nucleai  testing.  Consei]nenlly  catalytic 
destruction  ol  o/one  by  NO^  .eeneialed  dm ing  nuclear 
tests  could  nol  be  the  only  contiolhng  mechanisin 
alfecling  o/one  dining  till'  peiiod. 

2.  (),  and  l’(  .A's.  Theie  is  no  appaient  dose 
lesponse  lelationship  eilhei  m teinis  of  die  lelative 
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iiUciiMlies  of  the  PC'A's  1)1  I'-JOCi,  and 

or  In  lerins  ol  ilio  colatiliido  ol  llio  (),  dala. 

.V  r oirolalK)!)  ol  (),  with  snnspols.  Ono  nuMliod 
in  wliidi  o/ono  inighl  hf  atTcclcd  h\  solar  ai.livu) 
adopts  the  lollowini;  atjnnnont  (Rudcrnian  and 
( hanihoilain  An  ineiease  in  sunspot  ae'ivity 

iiKieast'i'  solai  wind  • decreased  eosniie  ta>  tinx 
into  the  ii[)per  atmosphere  — decreased  NO  piodiiction 
— incieased  o/one.  IIiis  theoiy  is  not  consistent  with 
the  ohserved  variations  ol  total  o/one  O,  oscillations 
decreased  m amplinidc  as  the  sunspot  aniphlude  in- 
creased lip  to  l'h)V)  lotal  incieased  lioiii  apeak 
sunspot  niaxiniiiin  ol  I '>5"’  to  the  vei-  low  siins(ioi 
inaMiiiiini  ol  l‘l(i’i  Heloie  l‘>54.  (),  oscillalioiis 
lajteed  sunspots  h\  ahimt  ihiee  ye.iis.  altei  I'h','  the 
variations  appeal  almost  in  [iha.e  O,  ainphliides 
exceed  ainphliides  coiiipnled  lioiii  ihe  Kiideiinaii  and 
( haniherlain  (14'’s)  iheoix  h\  ahv  nl  tenlold  There 
are  ralhei  lame  heinispheiic  dillerences  in  the  appaieni 
1 I -\  ear  iiscillaiions  in  O . 

4 Seasonal  chanees  in  lotal  (),  Hie  aiea  with 
the  hesi  lone-peiiod  (),  dala.  Western  I nrope.  shows 
between  H>55  and  HPO  a 14'.  inciease  in  winter  O, 
and  onis  4'  in  snnimei  O,.  This  sniiccsls  processes 
other  than  pholocheniisirv 

' Iheli  O,  niixnii;  ratios  in  polar  iiiulil.  Heath 
(H>^4)  repotted  III  \ observations  ol  a 4-nibar  O, 
mass  mixnia  lalio  ot  7.'  ppm  at  7()°N  on  2 .lannaiv 
l')71.  44'  above  Ihe  niaximnm  ot  lb  ppm  normallv 
observed  m the  tiopic.il  stialospliere.  I nless  the  hmh 
O,  measnremeiil  is  due  to  inslriimenl  dilTicnlties.  this 
IS  hiehly  siieeestive  ol  a iionphotochemical  O,  sonice 
noi  included  in  any  cm  rent  models. 

b.  Total  O,  lollovvmt;  Aitiini:  ernplioii.  Hinrorm 
itlobal  mixmti  ol  the  estimated  l-MT  slialospheric 
inieclioii  ot  Cl  is  computed  (Ryan  and  .Mnkheriee 
|b75l  to  cause  a decline  in  total  0■^  ol  approximately 


57t . II  the  bulk  ol  Ihe  Cl  liom  Apiine  lemained  m 
Ihe  sontherii  hemispheie  as  did  the  bulk  ol  Ihe  Atninti 
diisl.  even  lar,uei  chaiiiies  in  O,  should  have  occnired 
there.  C)^  obseivalions  show,  il  anylhinj;.  an  inciease 
in  total  O,  lollowme  the  einplion 

Due  to  the  spaiseness  and  possible  nnreliahilitv  ol 
some  ol  the  dala.  most  ol  these  dillerences  cannol 
be  resolved  at  Ihe  present  lime.  We  have  indicated, 
however,  several  aieas  m which  more  dala  ol  a reliable 
nainie  would  be  veiv  helpinl. 

T.flect  of  Sudden  Stratospheric  Warmiiiits 
on  Global  O/one 

To  sniinlale  and  lacihlale  rnrihei  invest iual ion  ot 
the  possibility  that  variation  m the  iilobal  o/one 
.iinonnl  is  caused  by  sudden  stialospheiic  vvarmmiis 
(SSW's).  observed  SSVW  vveie  tabulated  by  yeai.dale. 
and  measure  ol  miensily . It  was  noted  that  in  Ihe 
winters  havme  major  S.SWs  (lbb2-b.s  and  1 bb'^-tiS ), 
total  O,  not  only  rose  to  hiehei  levels  soonei  but  also 
lemained  above  the  iniiioi  SSW  winter  oT  |dbb-b7 
thionehoni  the  rolhvwinii  spring  and  snmmei.  The  data 
since  Ibs.s  indicate  that  a minor  SSW'  is  the  normal 
silnalioii  in  the  northerii  hemisphere  and  only  the 
iiiajoi  vvaimings  aie  abnormal.  One  can  then 
hypothesize  lhal  the  O.  trend  since  IdSS.dovvn  until 
Ibbd  and  lip  since  IbbT,  is  due  to  major  SSW's  alone. 
Hiis  would  also  explain  why  the  trend  in  the  southern 
hemispheie  appears  to  be  iiiiich  weaker  or  even 
leversed.  since  levver  and  weaker  SSW's  occur  in  the 
soniherii  hemisphere,  ( hristie's  ( |47.sj  ()^  analysis  tits 
this  hy  pothesis  very  well,  as  shown  in  Tig.  21.  where 
Ihe  years  with  major  SSW's  have  been  circled.  The  fit 
with  O,  trend  analyses  oT  others  is  not  so  good,  but 
this  ajipears  to  be  a problem  iilagning  all  cnireiil 
theories.  1 alei  analysis  suggests  lhal  lale-season  major 
S.SW's  cnlmmaling  aTler  the  end  oT  .laniiary  have  little 
ellect  on  mean  annual  O,  tor  the  following  year.  e.g.. 
the  majoi  SSW  of  l‘)  .ianiiary  through  b Tebriiary 
057.  .Admittedly,  some  of  these  arguments  are 
speculative,  hut  il  is  hoped  lhal  these  ideas  will 
slinmlale  fiirtlier  reseaich  and  exchange  of  ideas  which 
might  clarify  this  possible  climatic  mechanism. 


Fig.  21.  finic  scries  of  annual  mean  values  of  o/one 
deviation.  Circles  indicate  years  vvilli  major  sudden 
slralospheriv  warmings.  lAfler  C'lirislie  1973.) 


2.7  POSSIBLE  CAUSE  OF  THE  PRE  1970  INCREASE  IN  OZONE 


riic  upwaiii  iroiul  m lolal  (),  iii  llio  iioiilicm 
lK'ini>iilKMC  liom  ahiHil  I'L''’’  In  !'*’()  iiujihl  lu- 
alliiliiitahlo  to  a cliaiipc  ui  the  uciicial  >.iKiilalli'ii 
riiosc  [nopnsmi;  Midi  c\plaMalll>ll^  (op..  Kniiihvi 
01  al.  C'liil/oii  liavo  ioikIoiI  in  ho 

iinnsiiooitio  as  tn  ulial  oliaiipos  liaso  nsoiiiioJ.  Miuo 
uo  aio  lacking  iii  imdorslaiKliiig  n lliis  aioa  I iiipinoal 
dodiiollnii  sooiiis  In  mdioalo  that  a woakoiiiiig  nl  llio 
lladlos  oiioiilallnii  and  onnsoi|iioiitl\  nl  llio  oxoliaiigo 
nl'  air  hotwoon  tlio  linpnsplioio  and  liio  si lalnsplioio 
luiiild  load  tn  an  nioioaso  in  slialnsplioi  u (),.  1 Ins 
onnollisinn  rosts  sliniigly  nii  llio  npininii  pios.iilino 
piinr  In  llio  iimoiling  nl  llio  oataislio  o\^lo  llial 

llio  irnpnsplioio  uas  llio  niaini  sink  Ini  slialnsplioiK 
(),.  ('iinonl  llionis  is  lOinaikahK  onnsisloni  in  slinuiiig 
dial  ill  llio  o\lialin('ioal  Inwoi  slialnsplioio.  O, 
onnoonlialinns  aro  in  onooss  nl  llinso  I'lodiolod  h> 
plininolioinioal  ot|iiililii  iiiin.  \ sUmdnun  in 
linpnsplioiio-slialnsf'lioiio  oxoliango  wniild  appoai  in 
allnw  ninro  liino  Inr  O,  in  lliis  loginn  In  dooa\  In 
iho  plininolioinioal  oi|iiilihrnini  onnoonlialinns.  hni 
dial  on'ool  may  well  bo  nvoipnooorod  h\  ndioi  ol'l'oois. 

slnwdnwii  ill  linpnspliorio-slialnspliono  oxoliango 
also  rcdiioos  llio  upward  rale  nf  liaiislor  nl  N,0  in 
llio  2-X-In-.10-kni  lovci.  wlioro  il  is  nxidl/od  In  NO  h\ 
Ol'l)).  and  allnws  iimro  linio  nil  die  wax  up  Ini  dio 
\t()  In  bo  plinindissnoialod.  prndiioing  N,  radior  lliaii 
NO.  .\nd  sooniidly . sinoo  llio  lladlex  ooll  is  drixoii  h\ 
Iropioal  onnvoolinn  xxliioh  in  Iiirii  iiiaiiilaiiis  iipxxard 
pressure  nn  llio  irnpioal  irnpnpaiiso.  aiix  xxoakoniiig  oan 
bo  e.xpcolod  In  allnw  a Inworiiig  and  wanning  nf  dio 
Irnpioal  Irnpnpaiiso.  This  in  liirn  wniild  alinxx  a liiglior 
mixing  ralln  nl'  IliO  llirniigli  die  onid  nap  nl'  llio 
Irnpioal  iropnpaiiso. 

Tims  a wcakoniiig  nt'  llio  Hadley  oiroidalinn  willi 
rosiillani  slnwiiig  nl'  die  trnpnspliorio-slralnsphorio 
oxoliaimo  has  at  loasi  throe  oniisoi|iionoos.  all  nt'w  hioh 
appear  In  load  in  inoroasos  in  Iho  siralnsphorio  rosorvnir 
nt  (),: 

1.  Il  roduoos  die  rale  nl' Iraiisl'oi  nl' (),  Irnm  die 
sinrage  region  nt'  Iho  Inwor  siralnsphoro  In  nndorlying 
woll-imxod  Irnpnsphoro  and  O^-desirnyiiig  siirlaoc 
bniindary  layer  lie.,  il  rodiioes  the  Irnpnpaiiso  value 
nl  the  nl  the  l-l)  model).  This  shniild  rosnh  bnlli 
in  a rodiiolinn  in  the  Irnpnsphorio  onnoeiili alinii  and 
siirlaoe  dostriiolinii  rale  nf  O^  and  an  iiioroaso  in  lolal 
Oj  ill  die  storage  rosorvnir  nf  the  lower  siralnsphoro. 

2,  It  slows  the  upward  lliix  nf  N,0.  thus  allnwing 
ninro  lime  for  phnindissnoialinii  and  rodiioing  Iho 
fractinii  that  is  nxidi/ed  by  Ol'l))  In  form  NO.  It  also 
prosiimably  allnws  ninro  tmio  for  npoialinn  nf  Iho 
iinknnwii  Irnpnsphorio  sink,  thus  assuring  dial  an  oven 
smaller  fraoiinn  survives  in  pindiioo  V)  in  the 
siralnsphoro. 


V II  alinxxs  a Inwoiino  and  xx.nmmg  it  Iho 
linpkal  I inpnpanso.  w hloh  in  linn  allnw  s a liighoi  1 1 ,0 
mixing  lalin  In  pass  ihiniigh  die  lin|Hoal  linpnpaiiso 
onId  liap  and  In  IIKIoasO  llio  hlimidllx  nl  iho 
slialnsplioio  this  in  Inin  lies  up  .i  l.iigoi  liaolinnnt 
sliainxphoin  N()^  .n  IINO^.  lodiwing  iho  liaolinii 
..apablo  nl  jiJaAino  O,  ,.alaK  lioallx 

llioio  lomams  die  nagging  pinbloin  dial  die 
pio  l‘1  (I  iiKioaso-  III  lolal  (),  appeals  in  haxo  boon 
.1  plioiiniiioiinn  niiK  nl  llio  iinillioiii  lloililsplioio  while 
dio  ll.idlox  ooll  IS  nsiialK  logaidod  as  a global 
plionninonnii  I Ins  is  iml  nooossaiilx  an  insiiininunlablo 
pinbloin  Wo  know  dial  die  nnidioin  hoinisphoio 
(winlo:l  lladlex  ooll  Is  slinngoi  lli.iii  Iho  sniilhoin 
hoinisphoio  (siiininoi  ooll)  Ixxioo  as  slinng  aooniding 
In  Rood  and  \'look  ll'K'')).  and  alinnsl  Ion  linios  as 
sliniio  aooniding  In  Dulsoh  (1'I72)  and  dial  die 
linpioal  linpnpaiisos  nl  both  lioiiiisplioros  vaix  in 
iinisnn  nn  .in  amuial  oxolo.  i.o  . highosi  and  onidosi 
111  Jaiinarx  and  lobriiarx  Ipiishod  up  bx  the 
mil  ihorn-homisphoio  ll.idloy  ooll)  and  InxxosI  and 
wainiosi  111  .liilx  and  .\ngnsi  xxhon  iho  xxoak 
sniilhoin-homisphoio  lladlex  ooll  is  npoiatixo  (Sniilh 
l‘)P.').  rills  so'onis  In  indioalo  dial  the  ll.idlox 
oiroiilalinn  and  Us  onnirni  nl  die  linpioal  Irnpnpaiiso 
and  Irnpnshoiio-stiatnsphorio  oxoliango  is  sirnngK 
asxnimoliio.  duo  prosiimably  tn  bndi  die  disli ibiillnn 
nf  land  and  sea  and  the  oooontrioilx  nf  iho  oarlh's 
nrbil;  iho  fninior  aniplifios  laliliidinal  lompoialiiio 
giadioiits  In  sironglhon  the  nnilhorn-honiisphoro 
ll.idloy  ooll.  and  die  latter  oaiisos  the  roooixod  sniar 
dux  111  January  In  o.xoood  dial  nl  .Inly  In  noailx  T7r. 
again  sirongl honing  iho  nnrlhoin  hoinisphoio  ooll.  Unis, 
rcsiriolinn  In  the  nnrlhorn  honiisphoro  nf  the  pro- 1 ‘>’0 
iiioroaso  in  (),  boonmos  an  addilinnal  aigiiinoni 
sup|inrling  a ro'alinnship  bolxvo'on  die  iiioroaso  and  the 
lrn|msphorio-sltalnsphoi  10  oxoliango  rate  sinoo  the 
latter  and  Iho  nnoloar  tost  hypnihosis  aro  Iho  niilx  nnos 
prnpnsod  xxliioh  xvniild  lead  In  a honiisphorio  as 
nppnsod  tn  a ginbal  phonnmonnn. 

Finally,  il  is  nooossaiy  tn  oslahhsh  dial  llioro  xvas 
a prngrossivo  xvoakoning  of  the  Hadley  oiroulalinn  and 
a sinxvdnwn  in  Irnpnsplioiio-stralnspliorio  oxoliango  I'm 
apprnxiiiialoly  a dooado  piini  in  H»70-‘’l.  Hus  is 
alloslod  In  by  the  fnllnwing  ovidonoo. 

I.  Angoll  and  Knrshnvor  (|U74)  Inimd  a 
5-tn-7-nibar  dooado  iiioroaso  in  Irnpioal  Irnpnpaiiso 
prossino  xxliioh  ap|K'ared  In  have  begun  .ibniil  |U57. 
■Assnoialod  with  il  xvas  a warming  nl  die  linpioal 
linpnpaiiso  al  Snigapnio  and  (ian  nf  apprnximalolx 
I K dooado.  xvhloh  they  siiggosl  d m.Jil  bo  lolalod  in 
Iho  Iiioroaso  in  II, ()  mixing  ratio  m die  Inxxoi 
siralnsphoro  nbsorvoil  by  Maslonbrnnk  I HI"!)  over 
Washmginn.  IX'  . holxvoon  l‘H,4  and  I 'HO. 


2.  Moasiiiomonls  of  Jimyc  layer  stralosplieie 
aerr)sr)lsi  by  hi)lli  aiibome  eolleetnrs  and  reinole  sensing 
repoi  led  a niinimuni  and  virtual  disappearance  of  the 
layer  in  mid-l')7l  (lux  et  al.  IdVdl.  While  the  inaior 
thinning  ol  the  layer  Iruin  Id(i3  to  this  lime  has  been 
related  to  the  Ml.  .Agung  and  subsequent  volcanic 
eruptions,  the  I')’’ I levels  represented  a significant 
decline  over  the  pre-.Agung  values  reported  by  Junge 
and  I riend.  This  is  even  more  apparent  after 
application  to  the  pre-.Aguiig  data  of  the  5-io-6-fold 
correction  for  inipactor  collection  efficiency  which 
now  .qipears  appropriate  (Cadle  and  (liams  l'^)75). 
Since  gaseous  sulfur  and  particulates,  'd:e  all  other 
iiiateiial  not  episodically  injected  into  the  stratosphere 
(as  by  volcanoes,  eic.t,  must  enter  through  the  Hadley 
ciictilalion,  the  pie-,\guiig-lo- |d7 1 ihimiing  of  the 
Junge  layer  suggests  a weakening  vif  the  Hadley 
cnciilalioii, 

,v  Winstaiiley  (l'*7,')  reported  lor  the  Sahel  strip 


2,8  WORK  IN  PROGRESS 

Temperaliire  Feedhaek  in  a Stralospheric  Model 

I’eilurbations  to  the  stralospheric  composition 
alfeci  the  slialospheric  lem|>eralure  profile  via  the 
svilar  and  Uuigwave  ladialioii  balance.  Changes  in 
leni|ierature  affect  chemical  reaction  rales,  which  in 
turn  teed  back  on  slialospheric  composition.  The 
teiii(K'ialure  leedback  mechanism  has  long  been 
,recogni/ed  as  potentially  nnporlani  in  miligaliiig 
Iv’rtiirbafions  to  the  o/one  conceniralion,  but  it  has 
only  recently  been  included  in  sliatospheric  iiKivlels, 
Barnell  el  al.  ( H>75)  have  investigated  the  lenipeialure 
dependerl^e  ol  o/one  conceiilratioii  neat  the 
stralopause,  and  they  showed  that  a decrease  in 
leniperalure  lends  to  increase  the  o/one  concentralion. 

In  order  to  investigate  the  effect  of  leniperalure 
Icx’dback  on  the  conceiilratioii  ol  o/.one,  I Ll.'s 
oiie-dimensional  Iranspoil-kinelics  model  was  coupled 
to  a stratospheric  radiative  Iransfer  model.  The 
lem|ieralure  profile  is  delermined  using  a radiative 
equilibrium  calculation  with  convective  adiustment. 
The  model  includes  solar  absorption  and  longwave 
interaction  by  (),,  11,0,  and  CO,  along  with  solar 
absorption  by  NO,.  The  techniques  adopted  for 
Irealing  longwave  radiative  liaiisfer  are  the  same  as 
those  described  by  Ramaiialhaii  (l‘)74).  .A  band 
absorptance  forimilalion  is  used  to  treat  the  d.6-pin 
band  of  ().,  and  the  luiulamental  and  several  weak 
bands  of  CO,  and  its  isotopes  in  the  l.'s-pni  region. 
An  emissivily  lorntulation  is  used  to  treat  longwave 
radiative  Iransfer  by  11,0.  Solar  absorption  by  O,  is 
treated  using  the  empirical  lorniulalion  given  by 
I md/en  and  Will  (I'i7.0.  Band  absorptance 


of  Africa  and  India  a coiilimimg  decline  in  sumniei 
nionsooiial  rainfall  since  the  late  I‘t2()'s,  apart  from 
a temporary  reprieve  in  the  early  IdSO's.  The  lecent 
drought  was  the  result  of  the  cumulative  efiects  of 
a decade  of  below-average  monsoon  rainfall  with  less 
and  less  rain  each  year  through  Id70. 

Since  tropical  convection  cells  penelialing  the 
tropical  liopopause  are  likely  to  decrease  even  more 
than  summer  ITC/  (mtertiopical  convergence  /onel 
and  moiisoonal  rainfall  in  general.  Winstaiiley 's  data 
also  indicates  a progressive  weakening  (from  about 
l')57  to  l‘)70)  of  the  Hadley  circulation  and.  thus, 
of  iropospheric-slralospheric  exchange. 

The  above,  of  course,  does  not  consiilule  proof,  but 
does  appear  to  indicate  rather  strongly  that  a 
progressive  1 0-to- 1 3-yeardong  weakening  of  the  Hadley 
cell  and  of  tioposplieric-slratospheric  exchange  is 
responsible  for  the  progressive  norlhern-hemispheie 
increases  in  total  observed  up  through  1970. 


' 

lot niulations  are  also  adopted  for  treating  solar 
absorption  by  11,0  and  CO,.  The  empirical 

formulation  of  Tuther  (1976)  is  used  for  solar 
absorption  by  NO,.  Solar  absorption  by  O,  and  NO, 
are  treated  independently  because  absoiption  by  these 
s)X'cies  is  weak  in  the  region  where  their  absorptnui 
bands  overlap.  Solar  ladiation  scattered  from  the 
troposphere  is  included  assuming  an  albedo  of  O.T. 

A single  cloud  layer  is  included  at  b.5  km  with  42' i 
cloud  cover  as  suggested  by  Cess  (1974).  The  lapse 
rate  within  the  tropospheie  is  assumed  to  he  wet 
adiabatic  (6..''  k km),  and  the  temperature  at  the 
earth's  surface  is  specilied  to  be  2KK  k.  The  tropopaiise 
height  and  the  temperature  profile  within  the 
stiatosidiere  are  coni|niled  by  the  model.  The  CO, 
mixing  latio  is  assumed  to  be  .'20  ppm  by  volume. 

The  lenipeialure  |irofile  geneialeil  by  tlie  ladialive 
eqiiilibrinm  model  is  shown  in  I ig.  22.  This  prot'ile 
is  similar  to  teiiipeiaune  profiles  observed  at  tropical 
latitudes.  The  tempeiature  profile  from  the  (.  .S'. 
SlaiiJarJ  Alinosplicir  (I9(,2).  which  lepiesents  a 
midlatitude  average,  is  shown  for  compaiison  since  this 
jirofile  is  frec|ueiitly  assumed  foi  transport-kinetics 
calculations  which  neglect  tempeiature  leedback.  The 
uniK’itiirbed  tempeiature  piofile.  which  we  will 
hereafter  refer  to  as  the  "'ambient."  is  similai  in  sha|ie 
to  the  standard  alniosphere  but  is  a few  degrees 
warmer  above  20  km.  The  higher  lemperature  is  related 
to  the  model-generaled  o/one  coticentialion.  Theie  is 
excellent  agreenieni  when  the  imdlatitiide  o/oiie 
ptolile  lioiii  the  I.  S.  SuiiuhirJ  Al)iu>sphcn- 
Siipplfimvls  (19(,(,)  IS  used  to  generate  an  ambienl 
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Fig.  22.  Tempcralufc  profile  derived  by  the  model  for 
unperturbed  eonditions  eompared  with  the  L'.  S.  Standard 
Atmosphere  temperature  profile. 


tciiijx’rulure  profile.  Large-scale  dynamical  processes 
are  neglected  in  the  determination  of  the  temperature 
profile.  The  processes  have  a significant  effect  on  the 
temperature  profile  near  the  tropopause  at 
midlaliludes.  Including  these  effects  would  likely 
lessen  the  temperature  gradient  above  the  tropopause. 

The  model  was  applied  to  two  perturbations  to  lest 
the  significance  of  temperature  feedback  on  o/one 
concentration.  The  perturbations  considered  were 
stratospheric  injection  of  NO^  and  tluorocarbon 
production. 

Stratospheric  Injection  of  NO^.  fhe  effect  of 
stratospheric  iniections  of  NO^  b\  super.sonic 
transports  (SST's)  on  o/one  concentration  was 
computed  for  four  cases; 

1.  No  temperature  feedback;  .V.  SlamlarJ 
Aliiiosplwn’  (Ibfi2)  temperature  profile  assumed. 

2.  No  tem|x-rature  feedback; ambient  temperature 
profile  from  radiative  cc|uilibrium  model  assumed. 

.T  With  temperature  feedback. 

4.  Teni|X’raturc  feedback  without  including  solar 
absoiption  by  NO,. 

The  first  case  is  an  assumption  frei|uently  made  in 
transport-kinetics  calculations.  The  second  case  is  the 
reference  condition  lor  comparison  with  the 
teni|X'rature  feedback  calculation,  which  is  the  third 
case.  The  fourth  case  is  included  to  indicate  the 
significance  of  neglecting  solar  absorption  by  NO,. 
These  calculations  arc  not  intended  to  imply  the  effect 


of  a particular  licet  of  SST's  since  they  do  not  include 
the  injection  of  other  species  such  as  H,0  and  SO,. 
Once  the  significance  of  temperature  feedback  is 
determined  for  the  NO^  perturbation  alone  (whicli  can 
be  compared  with  other  calculations),  then  the  case 
of  temperature  feedback  with  multiple  |x-rturbations 
will  be  considered. 

The  reduction  in  the  o/one  column  caused  by  the 
stratospheric  injection  of  NO^  is  shown  in  f'ig.  2.1  for 
the  case  with  temperature  feedback.  Hie  NO^  injection 
rate  as  NO,  represents  a global  average  value  uniformly 
distributed  over  a l-km-thick  layer  centered  at  the 
indicated  altitude.  The  o/one  reduction  is  a steady 
state  value  obtained  from  a time-dependent  calculation 
as  it  approached  et|uilibrium.  Relatively  large  levels  of 
injection  of  NO^  are  required  to  reduce  the  o/one 
column  more  than  a few  percent.  The  o/one  reduction 
values  in  I'ig.  2.1  are  smaller  than  previously  reported 
(Chang  Id7.‘i)  by  about  a factor  of  2 due  to  changes 
in  the  chemistry,  especially  the  Oil  + 110,  reaction 
rate.  The  reaction  rate  presently  used  for  Oil  + HO, 
is  2 X 10'*  * cnrVs.  whereas  previously  we  used  the 
upper  limit  2 X lo  ’®  ern^/s  (llampson  and  (iarvin 
l‘)75). 

The  ratio  of  the  o/one  reduction  for  each  given  case 
to  the  o/one  reduction  with  temperature  feedback  is 


^ I : I 1 1 1,  t,J I 

1 2 4 6 8 10  20 

NO^  injection 
rate  - 10^  kg/yr  as  NO2 


Fig.  23.  Calculated  o/one  reduction  vs  injection  rate  for 
17'  and  20-km  injection  heights. 


slu>wn  in  hi”.  24  for  the  17-kiii  injection.  Iliere  are 
signiricant  differences  between  the  various  cases  for 
small  NO^  iniection  rates.  I'lie  two  cases  witlrout 
temperature  feedback  demonstrate  the  sensitivity  of 
the  calculation  to  the  assumed  temperature  profile. 
The  o/one  reduction  using  the  S.  SruiiJurJ 
Aliunsplnrc  (l‘)o2)  temperature  profile  is  as  much  as 
l.T;  greater  than  for  the  ambient  temperature  profile, 
the  difference  decreasing  with  increasing  NO^  injection 
rate. 

The  o/one  reduction  with  temperature  feedback  but 
without  solar  absorption  by  NO,  is  substantially  less 
than  foi  the  other  cases.  Neglecting  solar  absorption 


by  NO,,  whivh  is  important  for  the  pertuibeil 
stratosphere,  leads  to  an  overestimate  of  the 
stratospheric  tempeiature  change  and.  consecpieiitly . 
an  overestimate  of  the  temperature  feedback 
effect. 

The  lesults  for  the  20-km  injection  case  are  shown 
in  l ig.  2.^.  Although  the  lange  ol  \0^  injection  lates 
is  the  same  as  in  the  previous  figure,  the  lange  of  o/one 
column  reduction  is  much  grealei.  The  lesults  are 
similar  lo  the  1 7-km  injection  case,  but  the  difference 
between  the  various  cases  is  not  rpiite  as  large.  The 
o/one  reduction  with  temperatuie  feedback  is  about 
>'~r  less  than  that  with  fixed  ambient  temperature  for 
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Fig.  24.  C'aKiilak-d  r.itio  of  o/oni'  rcdiielion  tor  a givi'n  caw  lo  o/onc  r.-dmlion  with  feedback  vs  NO  injection  rate  for  a 
l7-kni  injection  height.  Calculalioii  assumes  the  NO^  injection  is  tiniformlv  dislrihiited  in  a 1 km-lhick  layer  centered  at  the 
injccdon  height. 
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Fig.  25.  C'akulaled  ratio  of  o/oiic  reduction  for  a given  case  to  o/one  reduction  with  feedhack  vs  \()^  injection  rate  for  a 
2(>-km  injection  height. 


small  \()^  iMjCclii'ti  rales,  and  llic  ero.ssovci  oeeiirs 
near  I4'V  o/otie  reduelioii. 

Hiioroearbon  Production.  The  iiatural  atmosphere 
for  these  ealeulalions  contained  IKi  with  an  assumed 
mixing  ratio  of  lo  '^  by  vohmie  at  the  earth  s surface 
and  a eonslani  CCI^  flux  at  the  surface.  The 
tluorocarboiis  (TCI,  and  (TyCU  were  iiiltodueed  at 
the  surface  starling  in  Id.sO  and  running  through  I 
based  on  world  nuorocatbon  production  data 
(McCarthy  l‘)^4).  The  computed  reduction  in  the 
ivone  cohinin  resulting  from  the  atmospheric  release 
of  fluorocarbons  during  this  period  ol  time  is  shown 
111  big.  2h. 

The  reduction  in  the  o/one  column  at  the  end  ol 
l')75  is  lor  the  case  with  temperature  feedback. 

The  o/oiie  reduction  for  the  fixcd-lenipetalure 


calculation  using  the  ambient  temperature  profile  is  a 
factor  of  I .()'!  greater  than  for  the  case  with 
Icniiieialure  feedback.  With  the  1.  S.  Slaihlan! 
Annosphen-  (1%2)  temperature  profile,  the  o/.one 
reduclioii  is  a factor  ot  I IS  greater.  .As  in  llic  case 
of  NO^  iniection.  the  results  are  sensitive  to  the 
assumed  lemperalure  profile  foi  the  fixed-temperature 
calculations,  and  lem|K’raIure  feedback  signilieanlly 
mitigates  the  o/one  reduction.  Teniisei attire  leedback 
reduces  the  iiiaxinium  peieentage  change  in  o/one. 
which  occurs  near  40  km.  The  maximum  change  in 
absolute  concentration  occurs  at  2H  km 
Summary  of  Temperature  Feedback  Investigations. 
Tlie  efiect  ol  teniperatuie  feedback  on  o/one 
concentration  has  been  investigated  for  two 
pel  tuibations  to  the  atmospheric  composition. 


rom|n.Taluie  Iccdback  was  I'oLmd  to  have  a similai 
etleel  lor  holh  the  sliatosplieiie  iiiieelioii  ol  N()^  and 
lliioroea] lion  piodiietioii.  In  boili  eases  leniperature 
leedhaek  bad  a 5-10' I icsloiini;  elTeet  on  the  cvone 
eolnmn  loi  small  o/one  leduetions.  Comparint; 
ealenlallons  based  on  eilber  rixed-temperalnie  |imfiles. 
teinpeiainie  leedhaek.  or  lemperatnre  leedhaek 
ne^leelin^  solai  absorption  by  NO,  revealed  these 
additional  important  points: 

1.  lixed-lemperaUire  ealenlalions  are  ipiite 

sensitive  to  the  assumed  tempeiatnre  prol'ile. 

2.  Calenlations  ol  .A()^  due  to  NO^  inieetion 
alone  inelndine  lemiieiatnre  leedhaek  may  be  as  mneh 
as  20' i less  than  those  using  tire  .S'.  Standard 
Amuisplurc  (I0(r2t  temperature  piolile. 

Cale'ilalions  that  negleet  solar  absorption  b\ 
NO,  sigm  .antly  overestimate  the  elTeet  ol 

temperature  leedhaek  lor  |>ertnrbations  involving  NO^ 
niieelion. 

As  the  ebemieai  kineties  and  tran.s|iorl  models  of 
the  stialosphere  eonlinne  to  improve  in  detail,  the 
elTeets  dne  to  temperature  leedhaek  will  beeonie 
ineieasmgly  important  and  mtist  be  eonsidered  in  any 
a.ssessment  ol  pr'ienlial  perinrbalions  to  the 
atmospherie  vi'inposition. 


Variation  of  Tropospherie  Minor  Species 

VVe  aie  examining  the  latitudinal  and  dininal 
vaiialions  of  tropospherie  minor  species  with  pai  lienlai 
atiention  to  Oil.  This  is  an  initial  elToil  lo  impioce 
the  transport-kmeties  model's  ea[iabilil>  lo  lie.it 
tropospheiie  residence  limes  loi  modeling  ol  aiieiall 
injeetions  and  also  loi  possible  replacement 
ehlorollnorometlianes  in  aeiosols  and  lel i igeiatioii . 
.Aeeoriling  to  model  ealenlalions.  most  ol  ihe 
tiopospherie  Oil  is  eonlained  within  the  tropical 
region  latiliide).  Oiiiinal  ealeiilatioiis  mdieate  a 

sinusoidal  variation  in  the  daylimo  eoneentiation  ol 
OH  at  tropospherie  allitiides.  willi  llie  ma.ximtim 
eoneentiation  oeemring  ai  noon.  Very  small  nighttime 
OH  eoneentiations  are  predicted. 

.Since  the  HO^  species  (OH.  HO,,  and  11,0,)  have 
very  slioit  lilelimes  in  the  liopospheie  (on  the  order 
of  a leu  seconds  lor  OH  and  10'^  seconds  lor  HO 

s 

as  calculated  by  the  model),  it  is  possible  lo  caleiilale 
the  alliliidinal  and  latiuidmal  varialion  of  OH  in  (he 
troposphere  willi  the  I -I)  model  uilhoiit  considering 
hori/ontal  transport.  I’lie  piimars  roimalion  of  OH  m 
the  tiopospbere  comes  throiigh  the  leactions 
0('D)  -t-  11,0  and  0('m  + Cll^.  while  the  loss  of 
OH  is  eonirolled  by  OH  + CO  and  OH  -t  CII  , I’rofiles 


1950  1955  1960  1965  1970  1975  1980 

Year 

Fig.  76.  (alailati’U  thungi’  in  o/oiic  cnliimn  duo  lo  t'luor(KarlK)n  prodiiclion  bolwcen  1950  and  1976. 
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Fig.  27.  Ditiriially  and  seasonally  averaged  values  ol  OH  eoneenlralion  vs  lalllude  and  height  in  the  troposphere.  Contour  values 
arc  concentrations  (!(/'  molecules  cm  ). 


ol  te'iiiiK'i .iluic.  11,0.  (lij.  CO.  anil  O,  (llmnigli 
pliololysts  to  inovidc  ()(M)|)  were  spcviiicii  in  llio 
nioJel  jtui  vailed  vvlili  lalllude  aeeoidine  lo  available 
nieasurenienls.  I lie  latiuidinal  disliibulion  ol  vvalei 
vapor  was  taken  lioir.  Ooil  and  l{asiiuissen  ( l')7l  ).  the 
(),  data  are  tioin  Dtilseh  (1‘Xi'M.  llie  lemperaune 
proliles  veisus  lalllude  are  lioin  I otiis  ( l''l"’4).  and  the 
CO  data  are  trom  Seller  (1474).  1 he  model  was  iheii 
used  lo  eakulate  the  diuinal  Oil  eoneenli ations  hu 
various  seasons  and  laliludes,  .A  eoin)iaiison  of  these 
ealeulations  wilh  the  available  ineasuied  Oil 
disitibulion  derived  Ironi  the  model  is  shown  in  I ig. 
21.  As  iiidiealed  above.  7|  ' o|  the  lro|iospheik  Oil 
IS  contained  wilhiii  the  band  .C()  S vO  N, 

The  derived  C>tl  coiieenlralions  were  used  lo 
delenmne  Ihe  liopospheiic  ehemical  hletime  loi  a 
niimbei  ol  niioioearboii  species  based  on  chemical 
degradation  via  Oil  ladical  The  calculaled  liletimes  aie 
shown  in  lable  10.  The  leaclion  tales  lot  ihese  species 
vvilh  OH  aie  based  on  mcasiiiemenis  by  Davis  ( l'>7(i). 
The  liopospheric  residence  lime  ol  these  species  is 


Table  10.  Estimated  tropospheric  chemical  lifetimes 

for  some  fluorocarbons  

rtuoroearbon  l.ife'liine'  (years) 


tii^n 

0.4 

CMjdj 

0.2 

( ( ijii 

0.2 

( M^Hr 

0.5 

( lljK  1 

0.5 

( Ml  1 

4.f> 

( III  ( Ij 

0.7 

( ll,(  1 j(  1 

7.1 

‘2^4 

0.2 

1.4 

impoilani  in  delei mining  vvhelhei  their  rmuie  use  will 
lead  to  an  im|\iet  on  stralosphetic  o/one. 
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Model  SciiMti\ity  I o Tropospheric  N,0  COiiceiUmlioiis 

We  luise  mvesliejled  the  elTeci  ol  doiil'liiii:  the 
Iroposphetic  N,C)  coiicentialiotis  iii  oui  siaiulaid 
model  and  m models  designed  lo  ma\imi/e  oi 
mimmi/e  the  ertcci  ol  an  N-,()  increase.  We  lonnd 
o/one  decreases  of  2.(V.  (minimum  reiinction).  7.5'' 
(standard),  and  I 1 .7''!  (maMimnn  reduction)  for  the 
Chani!  (l‘>74)  edriy  diflnsion  profile,  riuis.  the 
response  to  a doubling  id  liopospheric  N,()  was 
moderate,  and  the  ranue  of  uncertainties  was  not 
drastic. 

Nuclear  Weapons  Sensiti\it>  Studies 

riie  time-dependent  response  ol  the  model  to 
massive  injections  of  N()^  such  as  might  result  Ifom 
the  explosion  of  hundreds  lo  thousands  ol  megaloiis 
of  nuclear  weapons  was  calculated.  This  woih  is  still 
incomplete  and  will  not  be  described  in  detail  at  this 
time.  Ilowevei.  one  significant  lindingwas  that  a pulse 
injeclioii  of  NO^  in  ihe  lowci  slialospheie  coinpaiable 
to  that  expected  iVont  a lew  hundied  megatons  total 
yield  Using  l-to-.'-Ml  devices  lesulted  in  a tiansieni 
increase  in  total  o/one  when  models  .A|.  .At  ( |.  C,. 
C I , or  ('  j (described  in  Sec.  2-1  of  this  leporl)  were 
used.  This  suggests  that  single  test  clouds  might  result 
in  local  o/one  increases  after  a few  days,  malsnig  the 
results  of  such  studies  even  harder  to  interpicl  as  tests 
of  theory.  W'e  are  currently  regenerating  the  coding 
necessary  to  study  single  clouds. 

Two-Dimensional  Transport-Kinetics 
Model  Deselopinent 

During  the  last  year,  numerous  changes  weie  made 
111  the  tsso-dimensional  iranspoi l-kinelies  model  ol  the 
earth's  atmosphere,  in  the  mean  wind  and  eddy 
diffusion  coefficient  data,  in  the  chemical  idnctics 
data,  and  in  the  various  computer  codes  used  lo 
process  input  and  output  data.  Numerous  runs  have 
been  made  lo  vciify  the  behavior  ol  the  lianspoit- 
kinetics  model,  and  to  verify  the  lianspoit  model  by 
compaiison  with  experimental  data  for  several  passive 
tracers.  Cm  enl  woik  is  under  way  lo  improve  the 
accuracy  of  the  laiiludmal  deix'iidence  of  pholo- 
dissi  rciation. 

The  seasonal  average  mean  wiml  data  of  l.ouis  was 
fitted  with  I ourier  scries  to  provide  a lull  global 
linie-depcndcnl  mean  wind.  The  eddy  ilifliision  ihila 
ol  I uthci  was  mollified  lo  improve  accui.icy  in  the 
slialosphcre.  (ilobal  miilliplieis  were  added  to  allow 
modilicalton  of  the  mean  wind  and  eddy  dilfusion 
cocificienis.  A capability  of  siinulaline  paiticiilale 
scllhng.  using  the  .Slokes-Cumiiiigham  loimiila,  was 
added  to  allow  com|iarisoii  ol  the  model  with 
cxpeiimciilal  measurements  lor  v.iiiotis  n.ittnal  and 
artificial  tracers  for  which  pailiculale  settling  may  he 
sienilic.ini . 


The  liansport  ei|uations  were  replaced  with  a new 
mass-consistent  foimulalion.  which  is  much  beltei 
behaved  when  the  mean  wind  field  is  not 
mass-consislcnl.  as  when  hori/oiital  oi  veilical 
components  are  modified  by  the  global  nuilliplieis  A 
variable  input  paramctei  was  added  which  allows 
choice  of  any  weighted  combination  ol  fnst-oidei 
(upstieam)  or  second-ordei  (central)  spatial 

dilTeiencing  to  model  adveclion  by  the  mean  wind, 
lliis  allows  testing  and  partial  control  ol  the  ellecis 
ol  mimerical  dilfusion  and  dispersion  An  option  to 
keep  concentrations  fixed  at  any  selected  altitude  was 
added.  Ihe  tipper  bouiidaiy  condition  was  modified 
to  impiove  stability  and  avoid  negative  concentrations 
which  could  result  liom  ceitaiii  combinations  of  steep 
giadienis  and  wind  direction  These  changes  weie 
tested  111  a tracei  veision  of  the  2-1)  model  with  no 
chemical  kinetics  and  a single  passive  tracer,  with  a 
variety  of  luiis  simulating  peiiods  up  to  three  yeais 
and  with  diffeieni  inilul  and  botmdaiy  conditions. 
Stability  and  accuiacy  weie  gieatly  improved,  and 
coriect  steady-stale  solutions  were  obtained  lor  those 
cases  with  known  solutions.  The  lull  liansport-kinelics 
model  was  lun  lo  sinnilalc  a six-month  period,  with 
encouraging  lesulls 

Ihe  ti.icei  code  was  used  lo  simulate  the 
distribution  ol  pluloi,4tnn-2.'S  icleased  by  the  buimip 
of  a S\  Al’-'fA  nuclear  power  souice  in  the  stratosphere 
in  Apiil  l‘i()4.  which  was  described  by  lelegadas 
(|4(>h).  k'aiious  combinations  ol  the  particle  si/e  and 
density  were  tised  lo  model  settling.  It  appears  that 
sellling  significanlly  affected  Ihe  actual  disliibulion. 
and  could  not  be  modeled  i|uantiialively  with  a single 
panicle  si/e. 

Ihe  liacer  model  was  also  used  lo  simulate  the 
distribution  of  excess  carbon- 1 4 in  the  atmosphere 
released  by  atmospheric  nuclear  tests.  The  several 
sinnihilions  began  with  the  measured  distribution  in 
Maich-May  l‘>(i.f.  reported  by  Telegadas  (D)7|).  and 
simulated  up  to  three  years  using  various  imilliplicrs 
for  components  of  the  mean  wind  and  eddy  diffusion 
coefficients.  The  model-piedicled  distribution  of 
excess  carbon-14  for  .Inly  l‘)04  is  shown  in  I jg.  2.S. 
and  the  measured  distribution  is  shown  in  Tig.  2‘). 
(iood  i|iiantitativc  agreement  was  obtained  for  the 
inaxinuiin  mixing  lalio  and  foi  the  integrated  total 
load  in  the  whole  atmosphere,  the  troposphere.  Ihe 
slialospheie.  and  in  the  norlhein  hemisphere.  All  were 
within  IX'l  of  the  measiiied  values.  Agreement  with 
measured  conceit iralions  at  pailiciilai  locations  was 
i|uahlalively  good  in  geneial  and  ipiantilatively  good 
at  lower  latitudes.  The  mam  diffciences  were 
iKCUirence  of  Ihe  maximum  mixing  ratio  at  about  25 
km  instead  ol  20  km.  and  concenlialion  contours 
displaced  toward  higher  alliliides  in  the  pohii  regions. 


Nil.'  Iv^l  ovi'i.ill  a.uioomciU  «as  oliiaiiK'i,l  hy 

rmilliplyinj;  I oiiis's  iiK'an  \mikI  Jala  In  0.5  aiul 
I ullii'i  s Imii/niual  oiKI_\  Jillusion  aK'Hicii’iils  by  0.5. 
\ dilloronl  \cilkal  piolilo  dI  iIio  \nlkal  aKIy 
Jillusioii  C1VII1CH.MH  vuili  sMiallci  values  in  llio 
slialiisplicic  nuiilil  puibablv  iciluve  Ilic  altiliiilc  orioi 
III  llic  maviiiuiiii  iiiiMipj  lalio  aiul  polai  voiieciilialioii 
eoiiuiiiis, 

lliv  liacoi  iiiDilel  was  also  mod  lo  siiiuilale  die 
disirilniliDn  ol  itmesieii- 1 S5,  as  deseiibed  in  Slebbins 
(|0(OI,  lesiilline  Iumii  a lelease  ol  aboni  10  M(  1 al 
lb  kin,  lO^N  laliinde.  on  AiiyiiM  15,  lOo.S,  Se'veial 
runs  were  nunle  beyinnina  wuli  die  disiiibuiion 
deseiibeil  loi  Deeenibei  l‘>(i,s.  l aii  aiueenieni  nidi  die 
ineasuied  disliibulions  and  nilli  lotus's  siinulalioiis 
was  obtanierl  loi  die  liisi  eiylii  nionilis  ol  die 


siimilalioii.  I III  loiigei  linies,  results  diveijied.  piobably 
due  lo  die  laek  ol  a model  loi  lainoiil  ol  die 
liiiijislen- 1 S5  III  die  iioposplieie.  I lie  best  lesults  loi 
die  Hist  lew  iiiondis  weie  obtained  widi 
liiuis's  mean  wind  values  lerlueed  in  hall.  Hie 
verlieal  proliles  ajiaiii  indiealed  dial  die  veilieal 
eddy  dillusion  eoellieienls  should  be  redueed  in 
the  siraiospheie. 

Hie  lull  lianspoi l-kiiielies  code  is  piesenllv  beine 
lesterl  with  new  ehenileal  leaelion  rales  and  a lull 
annual  varialion  ol'  die  solai  /eiildi  aiiule  loi  eaeh 
lalilude.  Kapid  ehanees  in  eoneenlialions  al  die 
bormdars  ol  the  leiniiiiatoi  iieai  the  winlei  pole  iiiav 
iet|uiie  ntodiliealions  in  die  model  lo  avoid  die  use 
ol  piohibilnel)  small  lime  steps  01  lailuie  ol  die 
rhUeieniial  et|iialion  solvei  lo  eoiiveiee. 


TRACER  18.50  MONTHS  563.07  IflYS  OR  T • 4.865E*G7 

PLOT  12  USER  SPECIFIED  CONTOURS  FROM  6.JR4E+01  TO  7.9585+02 
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Fig.  28.  DisIrilMilions  ol'  ewess  rarbon  14  for  .Inly  IVfi4  as  laluilalerl  with  Ihe  liUMlimensioiial  miMli'l. 


37 


lliicerl:iiiitio>.  in  llu*  Odd  Nitrogen  of  tiu- 

SlnilospluTi- 

111  ail  olloil  to  iiiiprovc  oiii  iiiKlcislaiuliii^  ol  1.1111C111 
stialiisptk.'iii.'  iiiiiilels,  iiiJikiiiiii  lln.'ii  iiiii.ci laiiilK's  amJ 
liiiutatloiis,  a levicw  was  iiiidorlakoii  ol  llio  ikUI 
iiilroi;i.'ii  hikltict  ol  llic  slialosplicn.'  as  coiiipiiU'i.1  h> 


tlic  iiuhIoIs.  I’Iiis  K'vicw  lias  ievcak’i.1  ( I I an  iiiKoKaiii 
consensus  iliat  N,()  -t  (x'l))  is  iho  piincipal.  it  not 
llic  only,  soince  of  sliatospliciic  NO^^.  (2)  csliinatcs 
ol  llic  sliciiittli  ol  llic  N,()soiiicc  iniblislicd  since  I ‘<70 
\ai>  by  nioic  than  an  oidci  ol  inaitnilndc  (0.2  to  4 
X 10^  inolc'Ciilcs.  cm" ‘s );( .^ ) additional  N()^  sonices 


Fig.  29. 

nunilKT 


Mcasiircil  iliMiibiilmii  of  excess  carlMiii  14  for  liine  Augiist  l')b4.  Values  are  in  units  of  10^ 
of  samples  ihe  value  is  haseil  on  is  shown  in  parentheses. 


aloins  |Ht  gram  of  air; 


"V  v:*'  . 
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uliicli  liave  hccii  proposed  ineUidc  cosmie  rays, 
downward  traiisprnl  t'roin  the  tlicmiosphoie.  oxidalion 
of  NII^  (which  allematively  may  be  an  NO^  sink), 
and  upward  transporl  of  NO,  through  the  tropopause. 
as  well  as  sunnised  but  yet  to  be  identilicd  sources. 


Aerosols  and  the  Radiation  Budget: 

One-Wavelength  Versus  Full-Flu,\  Calculations 

We  have  compared  one-wavelength  solar  radiation 
calculations  with  I'ull-nux  calculations  (0.d85  to  2.5 
pm)  in  order  to  assess  the  limitations  and  uncertainty 
of  the  results  when  only  one  wavelength  in  the  visible 
is  used.  Comparisons  were  made  for  a standard 
atmosphere  with  tropospheric  and  stratospheric 
aerosols.  Calculations  at  0.55  pm  greatly  overestimate 
■Mie  scattering,  but  good  agreement  with  the  I'ull-llu-X 
calculatioti  has  been  obtained  by  scaling  the  aerosol 
rrptical  depth  to  the  llu.x-wcighted  aerosol  optical 
depth. 

(n  the  past,  assessments  of  the  solar  radiative  effecis 
of  aerosols  have  been  made  using  one-wavelength 
calculations  in  the  visible  spectrum.  One-wavelength 
calculations  have  the  advantage  that  they  require  much 
less  computer  time  than  the  full-nu.x  calculations. 
Their  disadvantage  is  that  they  are  less  accurate  than 
the  lull-nux  calculations.  If  a method  could  be  found 
for  correcting  the  one-wavelength  calculation,  then 
extensive  calculations  on  the  solar  effects  of 
stratospheric  aerosols  could  be  undertaken  at  a 
reasonable  cost. 

Two  corrections  must  be  made  to  the 
one-wavelength  calculation.  Otie  correction  applies  to 
the  fact  that  the  radiative  properties  of  the  aerosol 
at  the  given  wavelength  (usually  0.55  pm)  are  not 
representaivc  of  the  full-flux  radiative  properties.  This 
can  be  corrected  to  a large  degree  for  weakly  absorbing 
aerosols  by  scaling  the  one-wavelength  results  for  the 
changed  lluxes  by  the  ratio  of  the  nux-weigbted 
aerosol  optical  depth  to  the  optical  depth  at  0.55  pm. 
The  second  correction  applies  to  the  difference  in  the 
solar  flux  incident  on  the  aerosol.  In  the  case  of 
stratospheric  aerosols,  the  total  downward  flux 
incident  from  above  is  accurate  to  within  a few  percent 
of  the  I'ull-tlux  value.  However,  the  one-wavelength 
calculation  overestimates  the  flux  reaching  the  earth's 
surlace.  ('onser|uently , the  error  in  the  flux  ineident 
on  the  aerosols  from  below  depends  upon  the  surlace 
albedo.  In  the  case  ol  tropospheric  aerosols  the 
downward  flux  and  the  upwaid  flux  ineident  on  the 
aerosols  are  overestnnateil  by  a|iproximalely  the  same 
•actor,  so  the  error  is  neaily  independent  ol  surface 
.ilberlo  ami  can  be  corrected  ellectively  with  a constant 
scaling  factor 

The  change  in  the  iipwaid-scatteied  flux  at  the  top 


of  the  atmosphere  due  to  a stiatosphei ic  aeiosol  layer 
between  Ifi  and  22  km  is  shown  m fig.  .lO  as  a 
function  of  solar  /.enith  angle  and  surlace  albedo.  I he 
solid  lines  correspond  to  the  full-flux  calculation,  and 
the  broken  lines  correspond  to  the  scaled 
one-wavelength  calculation  at  0.55  pm.  The 
one-wavelength  results  have  been  scaled  by  a factor 
of  0.75.  which  is  the  ratio  of  the  flux-weighted  aerosol 
optical  depth  (computed  using  the  total  downward 
flux  at  22  km)  to  the  optical  depth  at  0.55  pm  The 
aerosol  optical  depth  is  0.05;  a /ero-order  logarithmic 
size  distribution  (Toon  and  Pollack  1076)  is  assumed; 
and  the  wavelength-dependent  real  and  imaginary 
indices  of  refraction  are  taken  from  I’almcr  and 
Williams  (1075)  for  757,  II, .SO^.  The  scaled 
one-wavelength  calculation  du|ilicates  all  of  the  major 
features  of  the  full-flux  calculation  (e.g.,  regions  ot' 
increase  or  decrease  in  solar  heating  and  the  local 
maximum).  There  is  very  good  agreement  for  small 
surface  albedo. 

The  change  in  the  net  downward  tlux  at  12  km  due 
to  the  stratospheric  aerosol  layer  is  shown  m fig.  .T|. 
A decrease  in  the  net  flux  indicates  a decrease  in  the 
solar  heating  of  the  troposphere.  The  calculations  show 
that  this  particular  aerosol  layer  can  lead  to  an  mciease 
in  the  solar  heating  foi  small  solar  zenith  angles  and 
large  surface  albedo. 

The  significant  variation  of  the  change  in  solar 
heating  with  solar  zenith  angle  and  surface  albedo 
indicates  that  global  mean  calculations  in  which  a solai 
zenith  angle  of  60°  and  a surface  albedo  of  0..5  are 
assumed  may  vary  significantly  from  diurnal  average 
calculations,  f uture  assessments  of  the  solar  radiative 
effects  of  aerosols  should  take  into  account  diurnal, 
latitudinal,  and  seasonal  variations  (as  well  as  variations 
in  surface  albedo). 

Climatic  Effect  of  Contrails 

A scries  of  experiments  to  test  the  climatic  effect 
of  increased  cirrus  clouds  due  to  jet  contrails  is 
cuirently  under  way  using  the  two-dimensional  climate 
model  ZAM2.  One  experiment  consists  of  increasing 
the  computed  cirrus  cloudiness  at  all  latitudes  by  a 
factor  of  1.2.  Gary  Hunt  at  the  British  Meteorological 
Office  has  (lerformcd  a similar  experiment,  and  we  will 
compare  our  results  with  his.  Results  will  also  be 
compared  with  those  from  the  two-dimensional 
COMI  SA  model,  which  uses  a scries  of  latitudinally 
arranged  onc-dimensional  radiative  cr|tiilibrium  models 
similar  to  those  of  Manabe  and  Mollcr  (1'>6I). 
Another  experiment  attempts  to  simulate  the  global 
climatic  effect  of  increased  cirrus  along  a flight 
corridor.  In  this  experiment,  the  cirrus  cloud  cover 
between  latitudes  .TO  and  50°N  will  be  incrementally 
increased  by  0.10. 
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Surface  albedo 


Contours  are  change  in  upward  flux  at  40  km 
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Fig.  30.  Cak'ulatcil  change  in  upward  flux  at  40-km  altitude  as  a tiinclion  of  siirfaL'c  albedo  and  solar  /onilh  angle.  St>lid  turves 
are  full  flux  calculation,  dashed  curves  are  one-wavelength  (0.55-^ni)  calculation. 


3.  FUTURE  WORK 


l uiuro  work  will  be  dirceteil  lowrin.!  luilhcr  ilala  beeome  available.  I uuiie  irerturbnlioii  M iidies  will 

renncmeiil  and  application  ol  a variety  ol  speciali/ed  einpbasi/e  nniltiple  pertnibations  (e.p..  NO^  and  IliO 

nunterieal  nurdelv  and  snpportnie  researeh.  with  baekeiound  (1^). 

Refineinents  whieb  are  planned  Tor  the  Hie  two-dimensional  model  will  undergo  I'urtlier 

one-dimensional  transport-kineties  model  meliide  development  Additional  tracer  studies  are  planned  to 

ineorporating  the  elTeets  ot  multiple  moleeulai  reline  the  pieserihed  tiaitspoit.  and  simulation  studies 

scattering  and  surface  retlcction  on  pliotodissociatioii  aie  planned  to  validate  model  cbemistiy. 
rates,  incorporating  a paraineteri/ation  of  the  water  (alculations  ol  the  latitudinal  and  seasonal 

vapor  budget  into  the  model,  and  reniiiiig  the  vaiiations  ol  the  solai  ellects  ot  peiturbed 

tropospheric  cliemistiy  so  as  to  include  the  full  sitalosplieiic  species  will  be  conducted  using 

methane  cycle  and  smog  chemistry.  Tlie  model  will  coiicentration  profiles  generated  by  the 
be  used  for  a number  of  sensitivity  studies  to  assess  Iw o-dmiensional  model.  Special  studies  will  be  made 

the  el'fects  ol  these  relinemenis  as  well  as  the  effects  using  the  two-dimensional  climate  model  ZA.M2  to 

of  new  reaction  rates  or  photochemical  pathways  as  assess  the  climatic  sensitivity  to  these  perturbations. 
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